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The Size-Variation of Resistivity for Mercury and Tin 


By E. R. ANDREW 
Royal Society Mond Laboratory, Cambridge 


Communicated by D. Shoenberg; MS. received 8th September 1948 


> ABSTRACT. Measurements have beer made of the low temperature resistivity of mercury 


wires down to 6 uw diameter and of rolled tin foils down to 3 wu thickness. A continuous. 


/ increase of resistivity with decrease of diameter and thickness was found, and it is concluded 


that this was due to a shortening of the mean free path of the conduction electrons by 


| inelastic collision with the boundary surfaces of the metal. ‘The theory of this effect for a 


wire is examined ; a rigorous theory for a foil already exists (Fuchs 1938). By comparing 
the experimental results with the theory, values of the product /p are found for both metals ; 
(/is the mean free path and p the resistivity for the metal in bulk). These values of /p are in 
satisfactory agreement with those obtained from Pippard’s measurements on the anomalous. 
skin effect for the two metals. 


SiN ERO DUCT TON 
N the course of attempts to extend the scope of the investigation on the 
] intermediate state of superconductors (Andrew 1948) to specimens of smaller 
dimensions, a number of mercury wires down to 6, diameter and rolled tin 
foils down to 3 » thickness were prepared: in this paper the results of resistance 


| measurements on these specimens in the normal state are described. It had been 


noticed in the previous work that the normal resistivity increased with decrease of 
diameter. J.J. Thomson (1901) first suggested that inelastic boundary collisions. 
of conduction electrons in a metal might cause a reduction in mean free path and. 
thus an increase in resistivity ; this effect would of course become more pronounced 
as the dimensions of the conductor are decreased. The effect was found with 
evaporated metallic films by Hamburger (1931), Lovell (1936) and Appleyard and 
Lovell (1937), and for wires of small diameter by Eucken and Forster (1934) and 
Riedel (1937). The lowest temperature used previously was 14°x., but at liquid 
helium temperatures the bulk resistivity of pure metals such as mercury and tin 
becomes so small that the mean free path is of order 10~? cm. and thus the size 
effects should be noticeable for wires and foils of this order of diameter and 
thickness. By fitting the experimental results to theoretical relations for the size- 
dependence, values are obtained for the mean free path of the conduction elec- 


trons in the bulk metal. 


§2. EXPERIMENTS WITH MERCURY 
2.1. Experimental Details 
The material used was Hilger 11,023, stated by the makers to contain | part in 


10° of silver as the only measurable impurity. The wires were made by filling 
glass capillaries of appropriate bore with liquid metal and freezing them, electrical 
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connections being made to reservoirs at each end of the capillaries. For diameters 


greater than 20, a pressure difference of one atmosphere sufficed to overcome ~ 


surface tension effects, and the specimens were those used for measurements on 
the intermediate state of superconducting mercury which have already been 
described (Andrew 1948). 

For diameters less than 20 additional pressure was applied over one of the 
reservoirs. A pressurizer designed by A. B. Pippard proved useful for this 


purpose after various modifications and additions had been made; a section of this 


pressurizer is shown in Figure 1. Pure mercury was held in the distrene cup F 


G5 EF Db 
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Figure 1. Figure 4. 


Section of pressurizer. (a) Method of mounting tin foils. 
(b) Section of potential clamp. 


inside the brass pressure chamber G. The capillary Z, 4-5 cm. in length, was 
sealed into a fine hole H in the top of the chamber with black wax, its lower end 
dipping into the mercury. This seal, and also the brass rim K, proved an effective 
seal up to the highest pressure applied (60 atmospheres). The upper end of the 
capillary was sealed with distrene cement into a hole in the bottom of another dis- 
trene cup A containing mercury. ‘This cup was carried by the hollow cylindrical 
ebonite spacer C, in which longitudinal slots had been cut to enable ready access 
of liquid helium to the capillary. Compressed air was applied through the fine 
bore German silver tube E, by which the pressurizer was suspended in the Dewar 
flask, ‘The two platinum wires M and N served as current and potential leads to 
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‘the mercury in the upper vessel. Only one independent electrical connection 
vicould conveniently be made to the lower mercury pool, and was made by the 
» platinum wire L, hard-soldered to the brass, current and potential leads being 
‘connected to the outside. The resistance of the mercury in the normal state was 
sufficiently high for the end effect in the brass to be negligible. The specimens 
were frozen from the top (low pressure) end first, so that the volume reduction on 
») freezing could be taken up by liquid under pressure, without the thread breaking. 
| Although the capillary of smallest bore tried, 2 diameter, was successfully 
filled and the thread frozen, 6 proved to be the smallest diameter which could be 


» diameter specimens were annealed for several hours at —60°c. and cooled slowly 
9\thereafter. A sensitive indication of strain is given by the breadth of the temper- 
_ ature transition into the superconducting state which, as the data in Table 1 show, 
became greater as the diameter was reduced below 12 p. 

The glass capillaries were drawn down from wider tubes and microscope 
#) examination showed that the cross-section maintained its circularity after drawing. 
# As the resistivity of the mercury at room temperature is accurately known, resis- 
# tance served to give the best estimate of wire diameter, provided of course that the 
“specimen cross-section was uniform. As a check, the inside diameter of the 
yi capillaries was also examined by immersing them in a liquid of exactly the same 
li refractive index under a cover glass; the correct liquid was found using the 
/ mineralogist’s Becke line test and box of graded liquids. The inside diameter 
could then be measured directly at any point along the length of the capillary with a 
§ travelling microscope or a calibrated eyepiece. 
Resistance was measured with a potentiometer, using the improvements of 
) Kapitza and Milner (1937). Since mercury becomes superconducting, measure- 
% ments below 4-17°xK. were made in a longitudinal field greater than the critical 
s value. ‘The temperature of the liquid helium bath in which the specimens were 
7 immersed was determined from the vapour pressure using the 1932 Leiden scale of 
i temperature. 


a0 
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2.2. Experimental Results 


The normal resistances of all the wires at 79° k. and at various temperatures in 
i the liquid helium range are given in Table 1, where they are expressed as a fraction 
j of the room temperature (liquid) value. The resistance ratios are plotted in 
' Figures 2 and 3. The full lines are explained in §4. The C series of specimens 
| (Table 1) are those of narrow diameter prepared in the pressurizer, and the B series 
4 are those of wider diameter described previously. At 1-05°K. the bulk resistivity 
| was almost entirely independent of temperature and therefore nearly equal to the 
/“‘residual resistance”, while at 4-2°K. it was many times larger and so mostly 
“ideal resistance”. The low value of resistance ratio, 2-0 x 10~°, for the speci- 
men B6 at 1:05°x. is an indication of the high purity of the material. 

Corrections for the magneto-resistance effect and for anisotropy have been 
made to the measured values before compiling Table 1. For low resistivities such 
as were obtained with fat specimens at the lowest temperatures, the magneto- 
resistance effect was quite appreciable, especially as the applied magnetic field 
| required to destroy superconductivity was several hundred gauss; by making 
measurements over a range of higher fields a fairly accurate extrapolation to zero 
| field could be made. : 

The resistivity of solid mercury is markedly anisotropic. ‘The ratio of the 
| resistivities parallel and perpendicular to the principal axis is 0-76 (Sckell 1930), 
7-2 
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and does not vary appreciably with temperature down to 82°K., the lowest temper- 
ature investigated by Sckell. At79°x., the ratio of the solid resistivity to the room 
temperature liquid value is 5-5 x 10-? for a direction parallel to the principal axis, 
7-2 x10 perpendicular, and 6:6 x 10~* for a random polycrystal. It follows: 
from the figures in Table 1 that the thinnest pressurized specimens are probably 
monocrystalline with the principal axis approximately perpendicular to the wire 
axis. In fact, for all but three of the fatter specimens (B1, B2 and the second 
freezing of B4), the ratio lies between 7-0 and 7:3 x 10-, suggesting a preference for 
the crystallites to be oriented with their axes normal to that of the capillary. The 
three exceptions showed lower ratios, and so presumably had their axes oriented in 
other directions. Any geometrical effect due to imperfect freezing, such as the 
formation of cavities, would have produced higher ratios, and hence appears to 
be absent. . 


log d (cm) log d (cm) 


Figure 2. Figure 3. 


The resistivity of mercury wires plotted logarithmically against wire diameter d. 
Figure 2, for 4:2° x., 3:0° kK. and 1°5°x. Figure 3, for 3:5° k. and 2-5° K. 
The full lines are the theoretical curves of Fuchs for diffuse scattering (see § 4). 


Since all but the three specimens mentioned have approximately the same 
crystallographic orientation, a correction should be made to bring these three into 
line. The correction will, however, be necessarily somewhat uncertain since no — 
data are available for the anisotropy of the ideal or the residual resistance of solid _ 
mercury at liquid helium temperatures, and moreover, the effect of anisotropy — 
on the extra resistance due to surface collisions is not known. ‘The mean value of — 
the ratio of the resistance at 79° kK. to that at room temperature for all but the three — 
exceptions is 7-2 x 10-*. For these exceptions the resistances at liquid helium — 
temperatures have been multiplied up in the ratio of 7-2 x 10-? to the 79°K. 
resistance ratio for each specimen. This correction assumes that the anisotropy 
of resistance is independent of temperature, and in the absence of other inform- 
ation, this is the simplest assumption that can be made. In any case only 3 of the” 
15 mercury specimens are affected. | 


The significance of the results is discussed together with those for the thin tin | 
foils in § 4. 
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§3> EXPERIMENTS WITH TIN 
3.1. The Specimens 


The material used was Johnson and Matthey 2,356, stated by the makers to 
contain less than 0-004°% impurity. The ratio of resistance at 3-8° K. to that at 18° c. 
for the thickest specimen, 1-8 x 10~4, appears to be the lowest recorded for tin and 
indicates that the material is of high purity. Foils were rolled out using an 
accurately ground and carefully cleaned rolling mill; a range of thickness. 
32,000 . was covered. ‘The specimens were cut in the fem of rectangular strips _ 
using a razor blade for the thinner ones and a guillotine for the thicker ones. It 
was necessary that the width of the strip should be much greater than its thickness 
in order that the foil should approximate to an infinite plate for the purpose of 
comparison of the experimental results with the predictions of the theory. For 
thicknesses of 80 y or less the foils were 1 mm. wide and 3-5 cm. long. For foils of 
thickness 0-1-1 mm. the width was increased to 3mm. The thickest “foil”, 
2mm. thick, was made only 4 mm. wide in order to keep the resistance large enough - 
for accurate measurement at liquid helium temperatures. Since the mean free 
path turned out to he of the order of 70, the size effect for this thickness was 
almost negligible, so that the insufficient ratio of width to thickness was not 
important. 

The method of mounting is illustrated in Figure 4(a), p. 78. The strips were 
laid flat on the paxolin board AB, and fastened down at the ends with small brass. 
clamps C and D to which current leads were attached. A few millimetres from 
the ends, small paxolin clamps pressed on the foil bearing the potential leads. A 
section of one of these clamps is shown in Figure 4(d). A groove was cut in the 
base of the clamp parallel to the width of the foil, and a loop of 36s.w.cG. copper 
wire was fastened along the groove so as to protrude slightly below the plane of the 
paxolin (dotted line). This protrusion was then ground away until the wire 
presented a plane surface flush with the paxolin. ‘The force used in tightening 
the clamp was thus distributed over the whole underside and was not gece to the 
potential wire alone. 

The coeficient of thermal expansion of the paxolin in the direction of the tin 
strip was about two thirds of the value for tin. For the thinnest foils the small 
relative contraction of the order 0-005 cm. was easily taken up by the slight lack 
of tautness in mounting ; the thicker and therefore more rigid foils dragged slightly 
under the clamps which had only been screwed up lightly to avoid biting into the 
metal, as was indicated by the fact that they bowed out slightly on warming up. 
‘The sharpness of the temperature transitions into the superconducting state and 
of the longitudinal magnetic field transitions from the superconducting state to the 
normal state for all thicknesses of foil, examples of which are given in the next _ 
paper, indicated the absence of appreciable strain. 

The thickness of the foils was known roughly from the room temperature 
resistance, and the mean thickness was found with an accuracy of 1°% after the 
experiment by weighing the portion between the potential points. 


3.2. Experimental Results 


Since the Debye temperature is much higher for tin than for mercury, resist- 
ance below 4-2°x. is almost entirely residual, the temperature-dependent part 
being only about 1%. In order to cover a wide range of bulk resistivities, as ‘was. 
possible with mercury in the liquid helium range alone, it was necessary to make 
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4) measurements also in the liquid hydrogen range (14—-20°k.). Since, however, at 


| 14°«. the resistivity is already about 16 times its value at 4-2°K., it is unfortunate 


al that the interesting range 4-14°K. was not readily accessible. 


The normal resistances of the foils at 3-8° k., 13-9° x. and 20-4°k. are given in 


s) Table 2, where they are expressed as a fraction of the resistance at 18°c.; these 


# resistance ratios are plotted in Figures 5 and6. The full lines are explained in § 4. 


Table 2. Resistance Ratios for Tin Foils 


| (4) 

/} Specimen | Thickness () = = 
(1) (2) (3) Temperatures (° K.) | 

79 20-4 NOL) 3°8 

E7 1,950 11-7 0-203 100-0 30-4 1-80 
E 10 915 12-0 0-210 98-7 SOSH 2: \) Sa 
E9 487 12-1 0-206 100-0 30°3 222) 
E6 255 12-6 0-207 —~ — 2:20 
E8 133 Ws? 0-205 102-1 S027 2°61 
| Bs} aged 12-3 0-205 — _- 3-02 
Ea 53°5 1220 0-209 103-1 30°8 2:64 
E 12 30:7 12-6 0-207 103-1 $3140) 3-10 
als ise 11:9 0-206 104-1 339 7/ 5333 
E 14 9-9 Bor 0-205 — — 6°75 
E17 7-46 1223} 0:203 115-0 42-9 10:8 


| (3) Resistivity at 18° c. (microhm cm.). 
» (4) Resistance expressed as a fraction of value at 291° kK. ( x 104 at 20-4°, 13-9° and 3-8° K.). 


Values of the room temperature resistivity are also shown in Table 2. It will 


_ be seen that there is a general increase in the resistivity at 18°c. from 11-7 x 10-® 
_ ohmcm. for the thickest specimen to 13-7 x 10-¢ohmem. forthe thinnest. Bridg- 
/ man (1925) has measured the resistivity of single crystals of tin in various directions 
at 20°c., and found the values 14-3 and 9-9 x 10-* ohmem. for directions parallel 
and perpendicular to the tetrad axis, so that the value for a random polycrystal is 


11:1 10-*ohmcm. It seems, therefore, that the process of rolling causes the 
crystal grains of the polycrystalline material to have a preferred orientation of the 


_ tetrad axis towards the direction of rolling. 


As the figures in Table 2 show, there was no significant difference from one foil 
to another in the ratios of resistivity at 79° K. and at 291°K.; the small scatter of 
values is probably due to variations in the temperature of measurement from the 
nominal 79° k., due to changes of oxygen content in the liquid air. Since the mean 
free path of the conduction electrons at 79°K. turns out to be less than 1 cm., 
no size effect is to be expected, and therefore the constancy of the resistance ratio 

"may be interpreted as a constancy of anisotropy of resistivity from room temper- 
ature down to 79°. 

There are some indications from the data in Table 2 that the anisotropy remains 
much the same down to 14°x.; if it is still constant down to 3-8° kK. (which is rather 
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uncertain since the resistance is residual), then by plotting resistance ratios rather 
than absolute resistance, the effects of variations in crystal orientation are elimin- 
ated. This procedure has been adopted in Figures 5 and 6. 


§4. DISCUSSION 


The experimental results show that the resistivity increases for both tin and 
mercury as the specimen size is decreased, but the possibility must be considered 
that the effects are due to strain or impurity. The mercury is unlikely to have been 
appreciably contaminated in handling, particularly in a systematic manner, but 
strain was indicated by the superconducting properties for wires of diameter less 
than 12. The increase in resistivity on the other hand was quite appreciable for 
diameters as large as 100u._ For tin the superconducting properties in dica ted no 


2.2 Id=l 


3-0 2-0 10 
log d (cm) 


Figure 5. Figure 6. 


The resistivity of tin foils plotted logarithmically against foil thickness. 
Figure 5, for 38° k. Figure 6, for 13-9° k. and 20-4° k. 
The full lines are the theoretical curve of Fuchs for diffuse scattering («=0). 


appreciable strain, so that only impurity need be seriously considered. While it is 
true that the thinner foils are likely to pick up more impurity in rolling than the 
thicker ones, the amount is unlikely to account for the factor of ten increase in 
resistivity found with the thinnest foils. 

A test that may be applied to the hypothesis that the increase in resistivity is due 
to strain or impurity is to see whether Matthiessen’s rule is obeyed. As the tem- 
perature approaches 0° k., the ideal resistance tends to zero roughly as the fourth 
power of the temperature (actually a 4-2 power law fits the mercury results best). 
By plotting resistance as a function of temperature, the temperature-independent 
resistance may be readily found from the experimental results of Table 1, and is in 
fact very little lower than the value at 1:05°x. If there is no effect due to electron 
collisions with the boundaries, then this temperature-independent resistance is the 
residual value due to strain’and impurity alone. Subtraction from the actual 
resistance at each temperature then gives the ideal resistance as a function of 
temperature, and if the observed increase of resistance were due to impurity or 
strain, this ideal resistance should be independent of specimen diameter. The 
curves obtained after this subtraction are shown in Figure 7 for four mercury 
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» specimens, and it is seen that there is a size variation of this apparent ideal resist- 
sojance ; the increased resistance cannot therefore be explained away as due to strain 
Myif Matthiessen’s rule is valid, although strain may make a contribution. For the 
1 (tin foils also, the application of Matthiessen’s rule to the data in T'able 2 for the 
foils investigated at the lowest three temperatures, leads to values of apparent ideal 
ol resistance which are greater for the thinner foils; the increased resistance cannot 
) therefore be explained away as due to impurity for tin either, although additional 


Himpurity may make a contribution. The theory of the increase of resistivity due 


98 lead to a size-variation of the apparent ideal resistance in the manner found experi- 
9} mentally. 

The variation of the resistivity of thin films due to inelastic scattering at the 
ul film boundaries has been the subject of theoretical investigation by Thomson (1901) 


a : 
Temperature (°k) 


Figure 7. The temperature dependence of the Figure 8. The theoretical variation of 


apparent “ideal” resistance of mercury for resistivity with size of conductor. 
four wire diameters: (1) 6m, (2) 10-54, /=mean free path in bulk metal, 
(3) 46, (4) 790. d=plate thickness or wire diameter, 
€=proportion of electrons scattered elasti- 
cally. 


plate, Fuchs (1938). 

—+—+—-— plate, Lovell (1936), «=0. 

— — — —wire, using Lovell’s assump- 
tions, «=O. 

cd ede PERI wire, Nordheim’s relation, 

«=0. 


who first suggested the possibility, by Lovell (1936) and by Fuchs (193 8). Making 
certain simplifying assumptions, Lovell obtained the following relation between 
the resistivity p’ of the thin film and that of the bulk material p: 


Pe CiaHiniae & bite St belies, (1) 


where d is the film thickness and / is the mean free path of the conduction electrons 
| inthe bulk metal. Relation (1) is plotted logarithmically in Figure 8. Fuchs made 
_ a more rigorous calculation and solved the problem for the general case when a 
proportion « of the electrons is scattered elastically at the surface (Lovell had 
assumed «=0). Curves of Fuchs’ solution for «=0, 3, 9/10 are also shown in 
Figure 8. 

For the tin foils these theoretical relations of Fuchs should be directly 
applicable, and the curve for « =0 (diffuse scattering) has been fitted as well as 
possible to the experimental data from Table 2. In the absence of other inform- 
ation the simplest value to assume for ¢ is zero. The process of fitting yields 
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values at each temperature of the bulk mean free path / and of the ratio of the bulk 
resistivity p to the room temperature value. The product of these two quantities 
should be independent of temperature if « is independent of temperature. The 
curves for the three temperatures using a constant product /p = 2-0 x 10-4 ohm cm? 
are shown as the full lines in Figures 5 and 6, and indicate a satisfactory 
agreement between theory and experiment. The scatter of points is most notice- 
able at 3-8°x. for the fatter specimens, but since the scale is logarithmic this — 
corresponds to only small differences in (residual) resistivity of the order usually 
found between various samples of the same material. If the bulk residual resist- 
ivity is isotropic, rather than having the same anisotropy as the ideal resistivity as 
is assumed by the method of plotting, the effect would be to cause the points for the 
thinnest foils to lie slightly above the theoretical curve in Figure 5. 

The accuracy of the product /p for tin is estimated to be about 25% from the 
scatter of the points. If the scattering coefficient « were considerably different 
from the zero value assumed, the product Jp deduced would be greater. Figure 8. 
shows that the theoretical curve for « =4 is roughly similar to that for «=0 but 1s. 
displaced by a factor 2:5. For e = 9/10 the factor is about 14. 

The rigorous calculation for a wire corresponding to Fuchs’ calculation for a 
film has not been carried out. Eucken and Forster (1934) have used a relation 
due to Nordheim: 

pip=1 edd) — wae ee eee (2) 
where d is the wire diameter, and the coefficient « has a maximum value 8/37 for 
diffuse reflection. The relation is also plotted in Figure 8. Its derivation does. 
not seem to have been published, but appears to be based on a simple geometrical 
argument involving the assumption that collisions with the boundary may be 
described by an appropriate relaxation time for the process as for lattice collisions ; 
the validity of this assumption seems doubtful. From the form of relation (2), if 
a is independent of temperature, the additional resistance should appear as an 
increased residual resistance, since /p is independent of temperature; this is not in 
agreement with the results for mercury illustrated in Figure 7, or the corresponding 
results for tin.» Thus it seems that Nordheim’s relation is inadequate. The 
solution for a wire using Lovell’s assumptions has been worked out, as is shown in 
Figure 8; an outline of the calculation is given in the Appendix. 

It will be seen from Figure 8 that the curves for a plate and for a wire, using 
Lovell’s assumptions, are very similar, being laterally displaced along the log d/l 
axis by a roughly constant amount which corresponds to a ratio of about 3 in d/l 
for the p’/p values dealt with here. Pending the availability of a rigorous solution 
of the problem for a wire, it therefore seems reasonable to assume that it would. 
yield curves similar to those of Fuchs a for plate, but laterally displaced by the same 
ratio. Although such a mathematical guess should be treated with caution, it is. 
unlikely to be in error by more than a factor of two. Comparison of the results. 
for the tin foils using Fuchs’ curves, which are directly applicable to this case, with 
some scanty data for tin wires suggests that the ratio 3 is quite reasonable. * 

The curve of Fuchs for « =0 has accordingly been fitted as well as possible to 
the experimental data for mercury from Table 1; the curves for the various. 
temperatures using a constant product /p =1-5 x 10-!! ohmcm? are shown as the 
full lines in Figures 2 and 3. In view of the uncertainty associated with the factor 


* The thinnest tin wires were made using the pressurizer described in § 2 suitably modified for- 
use in a furnace. 
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») 3 which has been applied, an accuracy of better than 50% in Jp cannot be assumed.. 
it Values of /p for mercury and tin have been calculated by Pippard from his. 
Ie 1,200-Mc/s. high frequency measurements (1947), using the theoretical treatment 
i of the anomalous skin effect by Reuter and Sondheimer (1948). For mercury he 
Tf finds the values of the product to be2-0.and3-9 x 10-1! ohm cm? for diffuse and elastic 
©) scattering respectively, while for tin the values are 1-0 and 1-5 x 10“ ohmem?,. 
» each with an estimated accuracy of 50%. The coefficient of elastic scattering ¢ in 
sf, the D.c. case is not necessarily the same as that for the high-frequency case, since 
/7 in the latter case only electrons moving in directions almost parallel to the surface 
19 contribute to the current, while in the D.c. case electrons moving in all directions. 
} contribute. Bearing this in mind, and also the limited accuracy of both deter-. 
minations of the product Jp, the agreement is satisfactory, and leaves little doubt. 
_ that the size-variation of D.c. resistivity is due to inelastic scattering of the conduc-. 
tion electrons at the metal surface; the correlation with the high-frequency 
) measurements suggests that a low value of « characterizes the scattering. 
Although no theoretical value of the product Jp is available for tin or mercury,, 
i the formulae for a free electron metal may be used to check that the order of 
) magnitude of the experimental values is correct. Simple considerations show 
1) that the resistivity of a free electron metal is given by p=mw/ne?l, where e and. 
» mare the electronic charge and mass, 7 is the effective number of free electrons per’ 
cm? and vw is the velocity of the electrons at the top of the Fermi distribution.. 
Quantum theory shows that v=(h/m)(3n/87)', where h is Planck’s constant.. 
= Hence the product Jp is given by 


Wee? Gian ik (Et mel re Ate. (3) 


and rough values of m may be found by substituting in (3) the values of Jp deduced. 
from the present experiments. The values of m thus obtained correspond to 0-43 
electrons per atom for tin and 0-57 electrons per atom for mercury. Such values. 
are quite reasonable, and imply that the experimentally determined values of /p are 
of the right order of magnitude. 
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APA EIN DIX 


The assumptions made by Lovell (1936) in calculating the resistivity of thin 
films will also be made here for a wire. ‘The diameter of the wire d is considered 
sufficiently less than the bulk mean free path / for all free paths to be taken as. 
starting from the wire surface. The collisions with the surface are considered 
inelastic, so that there is a uniform density of electrons in all inward directions 
from the surface. 

Consider one point on the wire surface. Path directions from this point will 
be designated by the polar coordinates 6, 6; # determines the azimuth in a plane 
normal to the wire axis through the point, and @ is the angle of elevation of the path 
above the plane. The fraction of paths from the point lying within the element of 
solid angle defined by 0, 6+d0 and 4, 6+d¢ is cos@d6d$/27. ‘The distance r 
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along any path to the next intersection with the wire boundary is d cos p sec 0. 
Again, following Lovell, it is assumed that for r</ the free path is r, while for r>1 
the free path is /. Hence the modified mean free path I’ is given by 


ie {" : | (ies 7 deena i Ios 08 | dd, 


0 0 cos—1(@ GOs P) 


‘where a=d/l. Hence for a<1 


Vi 


bao <1 4a— = {B(a)+a°B(a)}, 
pice Nibines 7 


ee =1+0~2{aE (=) - -B(a)t. 
pleas). 7 a Gta 


where E and B are complete elliptic integrals. 

"The case a>1 is considered in order to complete the curve although the assumption 
made above that all paths commence on the wire surface no longer holds. ‘The 
integrals have been computed graphically for the purpose of drawing the curve 
in Figure 8. 


and fora>1 
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Critical Field Measurements on Superconducting Tin Foils 
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ABSTRACT. The critical fields of superconducting tin foils have been determined from 
resistance measurements in a longitudinal magnetic field. Pomeranchuk’s relation, 
As=H (1+ »’/d), between the foil critical field Hg and the foil thickness d was verified; : 
and values of the length 0’ were found as a function of temperature ; at 0° kK. the value of 2 
‘was estimated to be 1:2 x 10~* cm., rising rapidly close to the transition temperature. It is 
noticed that \’ is approximately twice the penetration depth A for tin at all temperatures. 
Assuming the validity of Ginsburg’s interpretation of \’ as the sum of A and the surface 
‘energy parameter f, it is seen that A and £ are approximately equal. 
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§1. INTRODUCTION 


| Yt was found by Pontius (1937) that the critical field of thin superconducting 
if lead wires is greater than the value H, for the metal in bulk. A more complete 
. investigation of the size-variation of critical field for mercury in the form of 
| evaporated films has been made by Appleyard, Bristow, London and Misener 
_ (1939). The present paper describes the results of a similar investigation for tin 
' using thin rolled foils. Since, however, the thinnest foil was much thicker than 
the penetration depth A of a magnetic field, the information derived was necessarily 
7 less comprehensive than that for mercury obtained by Appleyard et al., who used 
films of thickness of the order A. Experiments were also carried out using the 
thin pressurized mercury wires described in the previous paper (Andrew 1949), 
but no useful results were obtained since the small changes of critical field were 
masked by effects due to strain. 
Pomeranchuk (1938) has shown that if there is a difference in surface energy 
NH ,2/87 between a plate in the normal and the superconducting states, the critical. 
field for the plate is given by 


dig LAG) eee re eS ee (1) 
where d is the plate thickness, provided that d>’.. Ginsburg (1945) has shown 
' that if the penetration laws of F. and H. London are assumed valid, and the: 
penetration depth is assumed independent of field strength, then A’=A+§, 
where 8 =87(«,, —«,)/H,? and «,, «, are the surface energies per unit area between 
the normal metal and an LTSUEMSE and the superconducting metal and an insulator 
respectively. 

The experiments described here show that the relation between H; and d has. 
in fact the form of equation (1); values of X’ as a function of temperature are 
deduced, and their interpretation is discussed. 


§2. EXPERIMENTAL DETAILS 
2.1. The Specimens 


The material used was Johnson and Matthey tin No. 2,356. The report of the 
maker’s spectrographic analysis indicated less than 0:004°% impurity. The foils 
were rolled out using an accurately ground and carefully cleaned rolling mill. 
The specimens were cut with a razor blade in the form of rectangular strips 1 mm.. 
wide and 3-5cm. long. The method of mounting is described in the previous 
paper (Andrew 1949). ‘The sharpness of the temperature transitions, and of the: 
longitudinal magnetic field transitions, for all foil thicknesses, indicated the: 
absence of appreciable strain, as did also the residual resistance measurements. 
given in the previous paper. 

The specimens were mounted vertically in the appendage of a Dewar flask at: 
the centre of a Helmholtz pair of coils whose field was uniform within 0-1% over 
the length of the foil. Resistance was measured by a potentiometer method, 
using the improvements of Kapitza and Milner (1937). A measuring current of 
1 ma. was used for foils 3-4 1 thickness, while for the thicker foils the measuring 
current was increased roughly in proportion to the thickness. ‘The magnetic 
field due to this measuring current could be neglected since it was found that close: 
to the transition temperature, where its effect should be greatest, a current three 
times as large caused an apparent reduction in critical field of only 0-1 gauss. ‘The: 
thickness of the foils was known roughly from the room temperature resistance, and. 
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the mean thickness was found with an accuracy of about 1% after the experiment 
by weighing the portion between the potential points. 

The temperatures of the specimens wcre deduced from the vapour pressure of 
the liquid helium in which they were immersed, using the 1932 Leiden scale of 
temperature. 


2.2. Special Precautions 


For the thinnest foil (3) the change in critical field that has to be measured 
‘ranges from about 1 gauss when H, is 3 gauss to about 8 gauss when H, is 240 gauss. 
While such a change may be readily observed, a reasonable accuracy in the deter- 
mination of X’ calls for critical field measurements with an accuracy of a few tenths 
-of a gauss, and hence involves a number of precautions. It was not necessary of 
-course to have an absolute accuracy of this order, and in fact one could work 
-entirely in terms of the currents passed through the Helmholtz coils. Since the 
slope of the threshold curve is about 150 gauss/°K. at the transition temperature, 
falling off to about 75 gauss/°xK. at 1:6°K., very careful temperature control was 
-essential in order to achieve the desired accuracy, and this was conveniently 
obtained using the automatic device described elsewhere (Andrew 1948b); the 
error in critical field measurement due to temperature drift was not more than 
‘0-1 gauss except at the lowest temperature, 1-66° K., where an error of 0-3 gauss 
was possible. 

It was important that the applied field should be accurately parallel to the plane 
-of the foil. The field coils were mounted so that they could be rotated through + 5° 
about any horizontal axis. The position of the plane of the foil was known within 
a few degrees beforehand. Longitudinal resistance transitions were measured 
for several coil directions a few degrees on either side of the position which had 
been guessed as correct. From the transition at the highest field strengths the 
-correct position of the coils could be found within $°. An error of +1° caused 
the transition to be displaced to field values about 0:2°% lower, and of + 2° to field 
values about 0:8°% lower. 

Since the transition temperatures 7, for various specimens are never quite the 
same, a correction for the differences is necessary. Such a correction involves 
-some assumption concerning the functional form of the threshold field curves. 
For instance. it may be assumed that H, is a function of (7,— 7), or of T/T,. The 
former implies that the threshold field curves are geometrically similar for different 
specimens, but subject to small shifts along the temperature axis, while the latter 
-assumes that the values of critical field become equal at 0°K. Near the transition 
temperature, however, all likely assumptions will lead to the same correction, while 
at lower temperatures H, is much larger and dH,/dT is smaller, so that the 
fractional correction is in any case small. For practical purposes the former 
assumption was therefore made, and in order to avoid as far as possible the 
application of a numerical correction to the experimental results, measurements 
were made at temperatures 7’ such that A7( = 7,, —T) took certain predetermined 
-values. The small remaining correction that had to be applied, due to the fact that 
the AT valucs were not quite equal to the predetermined values, amounted to only 
a few tenths of a gauss. 

The experimental procedure was therefore as follows: (a) the temperature 
transition was measured (i.e. resistance against temperature in zero field); 
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(b) the coil system was exactly aligned at AT=0-53°x.; (c) the specimen was 
allowed to warm up above the transition temperature; (d) systematic measure- 
ments were made of resistance as a function of field at various AT. 

For the measurements of the temperature transition, the earth’s vertical field 
was compensated, butit was uncompensated in the critical field determinations and, 
therefore, its value, 0-4 gauss, had to be added to the applied fields. The earth’s 
horizontal field (about 0-15 gauss) was not compensated in any of the measure- 
ments; its direction was almost in the plane of the specimen, so that no serious 
error is likely to result from its neglect. 


2.3. Determination of the Critical Fields 


‘The value of the applied longitudinal field which was taken ds the critical field, 
was the value which restored half the normal resistance. ‘This requires some 
justification since previous workers have sometimes taken the field which produces 
the first trace of resistance (e.g. Pontius 1937). It is well known that unstrained 
single crystals of tin give sharp temperature and longitudinal field transitions, 
whilst those for polycrystalline material are usually about 0-02°k. and 3 gauss 
broad respectively. Since strain is known to affect the temperature and field 
transitions (Sizoo, de Haas and Onnes 1925), it may be supposed that the crystal 
grains of which the polycrystalline material is composed are each under stress at 
low temperatures, due to the anisotropy of thermal expansion coefficient for tin 
and an inability to anneal-at such temperatures. If one assumes a Gaussian 
distribution of stresses for the crystal grains, then one has a Gaussian distribution 
of transition temperatures and critical fields, since the effects of stress, being small, 
may be taken as linear. This information alone, however, does not suffice to 
determine the shape of the resistance transitions as the geometrical arrangement of 
the crystal grains is also of importance. If, however, we first consider the fictitious 
one-dimensional case of a metallic line consisting of a random distribution of 
elements, some superconducting and some normal, then it is obvious that the 
fractional resistance will be proportional to the fractional length of normal metal. 
If the transition temperature and the critical field have a Gaussian distribution, 
then the fraction « of normal resistance restored either at a temperature Tin a 
temperature transition or at a field / in a field transition will have the form 


x 


= —| SiON ID (2) 


a) 


| where, for a temperature transition, x=(T'—T,)/,/2T, and for a longitudinal 


field transition x =(h —H,)/\/2H), (T, is the critical temperature, H, the critical 
field, and 7, and H, are the standard deviations of the two Gaussian distributions). 

Temperature transitions for the four foils, of thickness 3-20, are shown in 
Figure 1 and equation (2) has been fitted to them as well as possible (broken lines). 
It will be seen that the theoretical form for the fictitious one-dimensional model 
fits the results fairly well. Etching showed that the crystal grains varied in size 
from 0-1 to 1mm.: thus for all the foils there was generally not more than one 
single crystal extending through the thickness of the foil, but there were generally 
two or three grains across the width. _ For one foil (6:94 ) there were many large 
crystal grains, several of which extended over the whole width of the foil (1 mm.). 
It is satisfactory that for this foil equation (2) fits particularly well (Curve (2), 
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Figure 1), confirming that it is a good approximation to the one-dimensional model. 


The slight deviations from equation (2) for the more usual case in which there are — 


several grains across any width may be accounted for qualitatively by simple 
arguments. 
Since the crystal grains extend right through the foil, each affords entry and 


exit for lines of force in a magnetic field, so that approximately the longitudinal _ 


magnetic field transitions should be geometrically similar to the temperature 
transition. This is in fact the case, as is illustrated for the 4-11 » foil in Figure 2 
close to the transition temperature, and continues to be the case with roughly the 
same width (or slightly narrower) at lower temperatures. ‘The transitions were 
about 3 gauss wide for all specimens and temperatures. 

There will be a family of critical field-temperature curves corresponding to the 
various values of 7, and H, in the Gaussian distributions. Any member of the 
family near the mean will serve our purpose since, when plotted asa function of AT, 
they will all coincide. The half-resistance point on the resistance—field curves has 
been chosen to determine the critical field values because this critical field is very 
close to the mean of the Gaussian distribution, and also because the field value can 
be determined most precisely at this point since the curves are steepest there. 


§3. EXPERIMENTAL RESULTS AND DISCUSSION 


The foils used had thicknesses 3-43, 4-11, 6-94 and 19-41. In order to obtain, 
in addition, critical field values close to the bulk values, a wire of diameter 0-26 mm. 
was used rather than a thick foil. A foil of this thickness would have a demagneti- 
zation coefficient of 0-005, causing an apparent lowering of the critical field by a 
fraction of this order; the value of for the 19-4 foil, the thickest used for these 
measurements, was 4 x 10-4, the effect of which was negligible. For the wire of 
diameter 0:26mm. and length 5cm., the value of m was 1-3 x 10-4, and hence its 
effect also was negligible. The wire was cast in a thin-walled capillary and the 
glass cracked off; potential leads of pure tin were fixed on with a pinhole flame. 

The values of the foil critical fields H, are shown in Table 1, and were plotted 
against the reciprocal of the thickness d,* as is shown for three values of AT in 
Figure 3; the fact that the relations are linear is a confirmation of (1). The slope 
\'H,, the intercept H,, and their probable errors were calculated by the method of 
least squares. From these quantities \’ was obtained as a function of AT; the 
values are shown in lable 1 and are plotted in Figure 4, the vertical lines indicating 
the calculated probable errors. ‘The smooth curve that has been drawn through 
the points indicates that ’, like the penetration depth A, tends to infinity at T, and 
becomes independent of temperature at low temperatures. The calculated values 
of H, are also shown in Table 1, and are in good agreement with other published 
values (de Haas and Engelkes 1937, Daunt and Mendelssohn 1937, Andrew 1948). 

Unfortunately absolute values of A are not available for comparison with the 
present values of A’, but Pippard (1947) has obtained relative values (A —Ag.3.x) 
Using the relation Aj. ,/A?=1—(T/T,)*, which fits the data for mercury well 
(Daunt, Miller, Pippard and Shoenberg 1948), Pippard’s results lead to the value 
Xo.3° x =(0-6 + 0-1) x 10-°cm. The measurements of Laurmann and Shoenberg 
(see Shoenberg 1948) are in good agreement with this estimate. The absolute 
values of A thus deduced are shown in Figure 4 (broken line). If Ginsburg’s 


* In the case of the wire, the appropriate value is twice the reciprocal of the diameter, as a 
cylinder has twice the surface area/volume ratio of a plate. 
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interpretation of A’ is correct, the differences \’ — A give the values of f, and these 
also are shown in Figure 4 (chain line). Smoothed values of 4’, together with 
values of A and f, are given in Table 2. 

The interesting feature which emerges, and which is shown clearly by Figure 4, 
is the near equality of Aand £ at all temperatures. The same feature was found for 


| ie mercury also by Ginsburg (1945), who calculated values of A and 8 from the mercury 
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Figure i. Temperature transitions of four tin foils; Figure 2. Resistance curves in a longi- 
thickness (1) 3-43 p, (2) 6:94 p, (3) 4:11 p, (4) 19-4 we. tudinal field for the 4:11 » foil. Mea- 
The broken lines are given by equation (2). suring current 1 ma.; T,=3-728° k. 
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Figure 3. Critical field values Hg of the Figure 4. The values of \’ deduced from the experi- 


foils for three values of AT, plotted ments, plotted against 4T (=T,,—T), together 
against the reciprocal of the thick- with the values of the penetration depth A, and 
ness d. of B (defined as \’—A). 


film results of Appleyard et al. (1939). The result of this calculation is, however, 
at variance with that of Désirant and Shoenberg (1948 a), who combined the results 
of their absolute determination of A for mercury (obtained from susceptibility 
measurements on colloidal spheres and on thin wires) with the critical field 
measurements of Appleyard et al. Doubt has however been cast on the reliability 
of the thin wire results (Shoenberg 1948); the discrepancy may therefore prove 
PROC. PHYS. SOC. LXII, 2——A 8 
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unfounded. In any case the interpretation of \’ as A+ assumes, as Ginsburg 
pointed out, that A is independent of magnetic field strength; this is an assumption 
which may well be untrue. The correct interpretation of A’ therefore awaits a 
more complete theory. It is of interest to note, however, that if Ginsburg’s 
interpretation is adopted, the values of 8 which have been derived here for tin are 


Table 1. Values of H,, H, and 0’ 


Values of Hy (gauss) 


E1 E4 E5 E2 ( He 
194 694 4-11 3-43 eae 
3-724), 334. 3-728 oer? 


ANID T SS) 
0-010 
0-020 
0-040 
0-070 
0-150 
0-230 
0°530 
0-780 
1-180 
1-600 
2-080 


RRR RR bh UYU 
NVYreFtFtofFanuoeo 


* DN’ in 107° cm. 


Table 2. Values of X’, A and 8 for Tin 


) (smoothed) A B 
APCK) —“0-%cm.) (10-Fem.)_~—(10-F em.) 

0-02 8-0 3-0 5-0 
0:04 5:4: 22 ScD 
0:07 3°8 hoy Ma 
0-10 SiatS 1S 1-65 
0-20 2°40 1:18 1-22 
0-50 — 1-69 0-82 0-87 
0-80 1:46 0-69 0-77 
1-20 133 0-62 0-71 
1:60 25 0-60 0:65 
2°10 1-20 0:58 0-62 


of the same order of magnitude as the values of A obtained from experiments on tin 
cylinders in the intermediate state (Désirant and Shoenberg 1948b, Andrew 
1948 a), where A=87a,,./H,.? and «,, is the surface energy per unit area at the 
boundary between normal and superconducting domains of metal. This means 
that «,,. and («,, —«,) are of the same order, as one might expect. 


REFERENCES 
ANDREW, E. R., 1948 a, Proc. Roy. Soc. A, 194, 80 ; 1948 b, ¥. Sci. Instrum. (in the press) ; 
1949, Proc. Phys. Soc. A, 62, 77. 


APPLEYARD, E. 'T’. S., Bristow, J. R., Lonpon, H., and Misengr, A. D., 1939, Proc. Roy. 
Soc. .A,.172, 540. 


ea 


Critical Field Measurements on Superconducting Tin Foils 95 


Davunt, J. G., and MENDELSsonN, K., 1937, Proc. Roy. Soc. A, 160, 127. 

Daunt, J. G., MIttsr, A. R., Prpparp, A. B., and SHOENBERG, D., 1948, Phys. Rev. 74, 842. 

DeéstranT, M., and SHOENBERG, D., 1948 a, Proc. Phys. Soc., 60, 413 ; 1948 b, Proc. Roy. 
Soc. A, 194, 63. 

GINSBURG, V., 1945, J. Phys. U.S.S.R., 9, 305. 


) DE Haas, W. J., and ENGELKEs, A: D., 1937, Physica, 4, 325. 


Kapitza, P., and MILner, C. J., 1937, ¥. Sct. Instrum., 14, 165. 

Lonpon, F. and H., 1935, Proc. Roy. Soc. A, 149, 71. 

PipparD, A. B., 1947, Proc. Roy. Soc. A, 191, 370. 

POMERANCHUK, I., 1938, see Superconductivity, by D. Shoenberg (Cambridge : University 
Press), p. 98. 

Pontius, R. B., 1937, Phil. Mag., 24, 787. 


_ SHOENBERG, D., 1948, Amsterdam Conference on Metals. 
- S1z00, C. J., DE Haas, W. J., and Onngs, H. K., 1925, Leiden Comm., 180 ¢, d. 


On the Radiation from Transient Light Sources 


By F. BOOTH 
H. H. Wills Physical Laboratory, University of Bristol 


MS. received 2nd February 1948, in amended form 14th August 1948 


ABSTRACT. A theoretical treatment is given of the problem of analysing the radiation 
from transient light sources in the form of clouds of incandescent particles. Recently 
a method of correlating this radiation with various physical parameters has been suggested 
by E. F. Caldin. In the present paper a critical examination of this method is made, 
and exact formulae derived for special models. 

Expressions for the radiation emitted by a cloud of incandescent particles, all at the 
same temperature, are obtained for the following cases : (i) a uniform spherical distribution 
of equal, perfectly emitting and absorbing particles; (ii) a uniform distribution of particles 
which are not perfect absorbers, but in which (a) radiation is reflected specularly from 
the particle surface, and (6) the particles reflect light diffusely; (iii) a uniform distribution 
of particles in a thin spherical shell. 

It is shown that the condition for the radiation from a cloud of’ black body particles 
to be equal in intensity to that from a black body of the same size and shape is that 
R?ap)< M, where R is the radius of the cloud, a the particle radius, pp the particle density 
and M the total mass of the cloud. The cooling of clouds by radiation transfer is also 
examined theoretically. 

The bearing of the results on Caldin’s method is discussed, and more exact methods 
of analysing the cooling of this form of light source sketched. 


§1. INTRODUCTION 


HE following paper is a theoretical investigation of some problems suggested 

| by an analysis due to Caldin (1946) of the radiation from transient light 

sources. The light sources Caldin considers are all of the same general 

type. A metal is ignited either in oxygen, or in the presence of an oxidizing agent ; 

owing to the large amount of heat liberated, the temperature of the reactants 

becomes very high and part of the energy is radiated as light. All the sources 

show the same type of light intensity-time curve—a rapid rise to the maximum 
intensity, followed by a slower fall in intensity. 

A general theory of transient light sources of this kind would be very compli- 
cated since three distinct problems would have to be solved: (a) determination of 
the behaviour with time of the pressure, temperature and velocity of various parts 
of the system as the reaction proceeds; (b) some knowledge of the kinetics of the 
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oxidation would be required so that the proportion of metal oxidized at any time 
and the temperature of the products could be found ; (c) evaluation of the radiation 
reaching an external point from the source. 

For practical purposes all we need to know is the quantity and quality of the 
luminous radiation as found by (c). But although problem (c) is quite distinct’ 
from problems (a) and (6), these must first be solved to find the radiation emitted. 
Unfortunately the two problems cannot be solved independently. For example, 
the dynamical behaviour will depend upon the reaction kinetics, and the rate of 
oxidation upon the pressure and temperature of the material of the light source, 
that is, upon variables determined by (a). At the present stage an exact solution 
of the whole behaviour seems out of the question; drastic assumptions must be 
made about the first two problems (a) and (4) and attention concentrated upon the 
remaining problem (c). This is the procedure proposed by Caldin (1946a, b). 

In Caldin’s first paper (1946a) the cooling of a body by radiation from the 
surface is considered. It is assumed that initially, and during the subsequent 
cooling, the body is of uniform temperature. By equating the rate of energy loss 
from the surface to the rate of reduction of the internal energy, expressions are 
derived for the total intensity of radiation, brightness and temperature as a function 
of thetime. Inthe second paper (1946 b) the spectral distribution of the energy is 
found as a function of the time, and, by integration, the total amount of energy 
radiated at a given wavelength is obtained as a function of the initial temperature. 
All the expressions obtained hold only for a ‘“‘ Planckian radiator’’—that is, for a 
radiator whose emissivity is independent of the wavelength and temperature. 

Now consider the transient light sources. Suppose that at the maximum 
intensity of the light emission the chemical reaction is complete, and that the 
material is of uniform temperature and at rest. Then if we assume that the material 
behaves as a “‘ Planckian radiator” of uniform temperature, the theory can be used 
to analyse the decay of the light intensity. For example, from the total energy 
liberated at any particular wavelength, the initial temperature can be found. 
The validity of the method rests, however, on whether or not the radiating material 
can be treated as a uniform temperature ‘‘ Planckian radiator’. From obser- 
vations on the spectrum and temperature it seems probable that the light radiated 
after the maximum intensity is due to temperature radiation from incandescent 
solid particles—very probably oxide particles. If we assume that there is little 
radial motion, problem (c) is resolved into the determination of the radiation from 
a roughly spherical cloud of fine particles, each cooling by radiation. 'The temper- 
ature and density are not, in general, uniform. If the problem of cooling is solved 
by Caldin’s method two basic assumptions are implied: (1) the temperature 
throughout the cloud is initially uniform, and remains so as it cools; (2) the radi- 
ation from the cloud is the same as it would be from a solid surface consisting of the 
radiating material and coincident with the boundary between the outer layers of 
the cloud and surrounding space. 

The second assumption may be restated as two rather simpler assumptions : 
(2’) if each solid particle in the cloud is a perfect absorber and emitter then the 
radiation from the cloud is that of a solid black body at the same temperature ; 
(2”) if the surface emissivity of each particle is « then the radiation from the cloud 
is that of a solid black body at the same temperature multiplied by e. 

In the following account, we shall find the conditions which these assumptions. 
impose on the physical properties of the material in the light source. 
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Se eal VIN ACE LOIN TOE VAS Sila LOIN =(2) 
2.1. Radiation from a Cloud of black-body Particles (Assumption (2')) 
We first consider assumption (2’), assuming (1) to be valid. Let R be the 


} radius of a spherical cloud composed of spherical particles of radius a, density po, 
and at temperature 7; the density distribution is assumed to be spherically 
y symmetrical. If NM(r) is the number of particles per cm? in a small region round a 
point P, distant r from the centre of the cloud, then Ny, the total number of 

< particles in the cloud is given by 


“RK 
No=4a| M(r)r2dr=3M/Anpa®, ee (2.1) 
~ 0 


*» where is the total mass of solid material in the cloud. ‘To examine assumption 


(2’) we require the total amount of light radiated from a cloud of this form in unit 


| time. It is first necessary to derive the fundamental equation for the transfer of 
* radiation. 


Let J denote the total intensity of radiation of all wavelengths at a point P, in 
directions. ‘Then if dA is a small area element at P in a plane perpendicular to s, 
the amount of radiation which crosses dA in unit time and spreads into a small solid 
angle dw containing direction sis J(P,s)dwdA. Let P’ bea point distant ds from 
P along directions. In traversing PP’, part of the radiation is lost by absorption, 
but there is also a gain from the radiation emitted by the particles. Assuming that 


(dA)? and ds are large compared with a,* the number of particles within the cylin- 


der formed by dA and ds will be NdAds. 


Since each particle is a perfect absorber the amount of light absorbed in unit 


‘time by the particles from the light in direction dw will be I(P,s)za?N(r) dA ds da. 


The radiation emitted into dw per second will be 47aN(r)dAdsB(T) dw/4n, 
where B(T) is the total energy liberated by unit area of the surface of a black body 
in unit time at temperature 7. From radiation theory 


i ie oe ape eee” 1 ee (2.2) 
where o is Stefan’s constant. 
Hence 
— [I(P’,s) —1(P,s)] dA dw =7a®N(P) dA ds dw[(oT4/7) —I(P, s)], 
or (PE s)alids—na-N@)(ol|7)—J(P;s)|. «5. (253) 


Now dI/ds =01/0s + 01/cot, where c is the velocity of light. Since we assume the 
cloud at rest and T constant, 0//dt=0, and (2.3) can be written 


CUO meyer pete ete A TE (2.4) 
where T= maar | Nas. pce (2.5) 


Integrating the equation along direction s from P, (s=s, or T=7,) to Pg (Ss =5q oF 
T=T2), we get 

I(74, 8) exp 7, — 1 (72, $) XPT, =(0T*/m)(exp7,—expT;). ...-.- (2.6) 
Now apply this equation to determine the radiation from the cloud. We first take 
‘N(P) uniform, but this restriction will be removed later. Thus 7=«s where 
a=na2N,. Let PQ denote a small cylindrical volume element of the sphere, 


* The only essential limitation this condition imposes, is thata<R. This is a fairly obvious 
restriction, for it would seem absurd to analyse the radiation transfer in the manner we intend, 
if all the material in the flash were concentrated in a few pieces. It transpires, however, that 
even in this case, the results of the method are correct. 
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whose axis is parallel to s and whose normal cross section is dA (Figure 1). The 
total light emitted from the end surface at Q in dw about direction s is [(Q, s)dA dw. 
From (2.6), 
I(P,s) —I(Q, s) exp (aPQ) =(oT*/a){1 —exp(«PQ)}. ...... (2.7) 
Clearly J(P,s) =0 since we assume no light reaches the sphere from outside, and 
so 
1(Q,s) =(¢T4/7)i1—exp(—ePQ)}, eee (2.8) 
If z is the distance of the axis of PQ from O, the centre of the sphere, dA =z dzd@ 
where 6 is the angle between the normal to PQ from O, and a reference line in the 
plane through O perpendicular tos. Hence from (2.8) the total radiation emitted 
into dw is 
ben Oct 
H =do(cT*/7) | | 2(1—exp {—2a4/(R*—2*)})dzd0, —...... (2.9) 
0-0 
which integrates to 
1 1 i 
4 72 sae cat <n ——_ 
dwaT*R E 7Ro2 + exp (—2Ra) (ze + sR) |- ee (2.10) 


Integrating over dw, the total radiation emitted becomes 


1 1 1 
aa 2 4 oe — ZR ae att 
H =47rR?*oT E 7p te (ze + RE) | 


= HEREC Ro), ye ae ey ees (2.11) 
where H*=47 Roofs eee ee (2.12) 
that is, the radiation from a solid black body sphere of radius R, and E(x) is given by 
% 1 i ee ae 
E(x) =1—- aa t & ( + 53) Neate (2.25) 


Clearly we may call E(x) the ‘‘ effective emissivity’ of the cloud. 

Consider the behaviour of E(x) as x varies: (i) as x— 00, E(x)—1; (ii) as 
x0, E(x) 4x/3; (i) dE(x)/dx =x-3[1 — e-?*(2x2+2x%+1)]. Since x is always 
positive e??>1+2x+2x? and so dE(x)/dx>0. Hence E(x) is always positive and 
increases steadily from 4x/3 to 1 as x moves from small to large values. Therefore 
the condition that H=H* or that assumption (2’) is fulfilled, is that Ra >1, or 
from (2.1), 


Rap, <9M/l6n. ee Eee (2.14) 

On the other hand, for Ra small or when 
Rapp > OM iba. © ee eee (2.15) 
H=4ReH*)3=4ra-Niol tee (2.16) 


The radiation is then simply the sum of the radiations emitted by each particle 
and none is absorbed. 

In subsequent work we shall frequently refer to (2.14) and (2.15) as E and E’ 
respectively. For brevity, when a cloud satisfies inequality E, we shall say it is an 
“KE” state, when inequality E’, in an ‘“‘E’”’ state. 

Now let us consider how the results of the analysis are modified if we remove 
the restriction that N(r) is constant. Equation (2.9) then becomes 


H=4oT | : | ~a(l —exp{—7(PQ)})dzd0, ...... (2.17) 
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‘) where 


(PO) =70" fe NGVAROy ee (2.18) 


To evaluate H accurately, the exact form of N(r) is required. But it seems very 
likely that condition (2.14) still applies for assumption (2’) to be valid, for it is 
easily derived for the case of a thin cloud in the form of a spherical shell, which is a 
probable density distribution for practical cases. Suppose the material is 
_ concentrated in a thin spherical shell of thickness A. Then H may be written 


: ron [- -R-A R 
H=4oT*| | 2(1 —exp {—7(PQ)}} t+ | z(1—exp {—7(PQ)})dz | da. 
4 0 0 ~ R-A 


Since A <R, the second term is small compared with the first and we have very 


nearly 
-R-A 
0 


H=4oT!| |“ 2( -exp {-2(PQ)}) dz aa. (D0) 
0 


For H = H* then, we must have, for all possible values of PQ, exp {—7(PQ)}<1 
or 


cht 1) Fs ch iD Aalto (2.21) 


Now in the range of integration 0 to R—A, r( PQ) has its minimum value when 
PQ lies along the diameter in the direction s. Hence (2.21) is satisfied if 


rem N Ge ee Nem ines (2.22) 


J R-A 
Since 47R? itm AN ar =No, (2.22) gives R®ap) <3M/87 which is equivalent to E. 


When 7(PQ) <1, it is easy to show that we can derive expression (2.16) for the 
total radiation emitted for in this case 


rR 2x -Q 
H=4oT tna? | | | N(r)zdz d d(PQ) =47a’oT*N,. 
0-0 P 


Summarizing these results, we see that for a spherical shell of radiating particles 
we may write 

eS ROE ee ee ee ae (2.23) 
where 


T) =7a" LS CN atl lat ti ee 2.24 
, R-A 


and E(x) ~ E(x) for both large and small x. 
Hence the main conclusions reached for the case of a uniform distribution 
still hold for a spherical shell. 


2.2. Radiation from a Cloud in which the Particles are not Perfect Absorbers or 
Perfect Emitters (Assumption (2")) 


We shall now examine assumption (2”). Suppose the absorption coefficient 
of the surface of the particle is independent of the wavelength. Then, by Kirch- 
hoff’s law, the absorption coefficient equals the emissivity «, where « is defined as 
the ratio of the total energy radiated by unit area of the surface in unit time to the 
energy radiated by a unit area of black body surface in unit time. To evaluate the 
radiation from the cloud, it is now necessary to take into account scattering of 
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radiation by the particles. For, consider a single particle, in the path of a parallel 
beam. If J is the intensity of the beam in direction s in dw, the total amount of 
light received from the beam in unit time is clearly Iva*dw. But the radiation 
absorbed in unit time is eJ7a2dw. Hence the amount scattered out of the beam by © 

the particle is (1—«)/7a?dw. Let a fraction 


»(P, 8’, 8) des’ 4a , ae here (2.25) 


of the total light scattered by a particle be scattered into directions within solid 
angle dw’ about the direction s’._ It follows that 


f(P,8ys)do' a4, = (| (2.26) 


where integration is over all directions s’. The more general form of equation (2.3), 
in which account is taken of scattering is 


oe natn | oT" HM ar SoS =i I(P, s)’7(P, s’ ,s) do |, 


The first term of the right-hand side represents radiation gained by emission, 
the second, radiation lost by absorption and scattering, and the third, the contri- 
bution due to radiation scattered zmto direction s. ‘The integration in the third 
term is again over all directions s’. 

If we neglect the light scattered into dw’ eg (2.27) becomes 


ae natn | SP _I(P, s)|, Bee (2.28) 


which is equivalent to (2.3) with o replaced by eo. Hence if (2.28) were satisfied, 
together with the condition E, assumption (2”) would be correct. But it is clearly 
illegitimate to reject the light scattered znto dw and retain the light scattered out in 
the derivation of (2.27). Hence (2”) is incorrect and gives too small a value for the 
emitted radiation. 

On the other hand, even if E is satisfied the radiation is Jess than it would be 
from a solid sphere emitting black-body radiation. For, suppose the cloud is 
placed in a black-body enclosure whose walls are at temperature 7. Let H,, be 
the radiation emitted by the material of the cloud which traverses the outer surface 
inunittime. ‘Then the total radiation received by the enclosure from the particles 
is H, =H,,+H, where H, represents the part of the energy from the walls of the 
enclosure which is scattered back to the walls by the cloud. But H, represents the 
energy emitted by a cloud of black-body material since in this case H,=0. Thus 
fore <1, H,,<H*. Hence it follows that the total energy radiated per second is 


H=4nR%'oT', (2.29) 


eee eee 


where «<e’<1. ‘To determine <’, which we may again call the “effective 
emissivity’? it is first necessary to determine y(P,s’,s). . This is, strictly, a matter 
for experiment. But y(P,s’,s) can be calculated for simple cases. For example, 
if the light scattered is reflected specularly from the particle surface 


y(P, s,s) =1. 
For diffuse reflection 


y(P, s’,s)=(8/37)(sind—fdcosd), esa (2.31) 
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| where ¢ is the angle between s ands’. Another law which might be considered is 
that for “ Rayleigh scattering’’, which is given by 


we, SS) = + ces S)i4., sty | he. (2132) 


But for this law to hold, a would have to be small compared with the light wave- 
9) length and it is unlikely that the particles in the light sources we are considering 
are so small. The evaluation of the expressions (2.30) and (2.31) is given in 
Appendix I. 

It now remains to calculate «<’ by substituting expressions (2.30) or (2.31) in 
equation (2.27) and solving for. Unfortunately, owing to the last term of equation 
{| (2.27), the total radiation from the sphere cannot be evaluated in the same manner 
& as in the case of (2.11). We can, however, give an approximate calculation, the 
1 results of which agree quite well with (2.11) in the case of « =1. 

i Let the angle between the radius vector drawn to any point in the cloud and 
? the directions be #0. Then 


ol ol sin @al 


ee Apo yt one aieiel cers (2.33) 
If we again assume JN is constant and write 
7=na"N(r), eee fe ee (2.34) 


then equation (2.27) becomes 
Clacin telco ES, (le) ‘ : aS 
Oe a | I(r, 0’)y(r, 0", 0)dw’ 1. ee ec (2.35) 


Let.J denote the total intensity of radiation and F the total flux of RIEL TOD | in the 
‘direction of increasing 7, then 


COs 


ONS | Ide =2n| I(r,8)sin9d8, saa (2.36) 
J0 
ee | Zoos een iC 6)sin@cos@d8. sss (2.37) 
0 
- We also define a quantity K by the relation 
ene i Feo 6 dw = 20 fe I(7,8)sin 8 cos?0d0. ss... (2.38) 
0 
Multiplying equation (2.35) by dw and integrating over all directions, we have 
dF ol ee etal : j : 
——- 0 6’, 0)dw' —J....... 2239 
ie 7. ~ | $sin bdo = 4eoT# + | deo | I(r, 6’)y(z, 8", 8) doo (2.39) 
Using (2.26) 
| de [4c 0’ )y(7, 0’, 8) dw’ = | T(z, 0’) dos’ | y(t, 0',0)duw=4rJ.  ...... (2.40) 
Also 


* 2" sind cos dd | So 
0 


ol 
Pd? [5 »9d0—2 [7 29 
J agsin dw =27 ap sin a7} Jsin ' 


Hence (2.39) becomes 
dFldr+2F/r=c(4oT*—J). seen (2.42) 
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Multiplying oa (2:55) phe cos@dw and integrating again, we get 


dK 1vol. (=) ; 
= 6’ 6’, 0) dw' —F 
ae are ag sin 6 cos 0 dw = a | cos 6 dw | I(r, 8’)y(z, 0, 8) dao 
Lok ee (2.43) 
Reversing the order of integration, we find 
| cos 6 dw { I(r, 8")y(r, 6, 6) des’ = | I(r, 6°) deo! | Ar", 0) cos6dw. 2... (2.44) 


Let ¢ denote the angle between s and s’, and ¢ the angle between the plane con- 
taining directions s’ and s and the plane containing direction s’ and the radius 
vector. ‘Then we can put dw=sin¢dd¢dy, and since 


cos@=cos¢cos@ +sin ¢ sin @’ cos ¢, 
we have 
rm 727 i. 
[ve 0’, 0) cos 6 dw =| | y (¢)(cos ¢ cos 6’ + sin &’ sin¢ cos #) sin ¢ dd dys 
0 0 


= —4nX cos 0’, me ee 


where X isaconstant. (¢) is given by one of the relations (2.30), (2.31) or (2.32). 


Integrating (2.45) we find X =0 for regular reflection or Rayleigh scattering, 
and X =4/9 for diffuse reflection. 
We also have 


ge sin 6 cos 6 dw =2z7 es sin? 6 cos @ dé 


ld 900 
=27 E sin? 6 cos 0 


* — ["1(3 sin @ cos? @ sin 8) ao | J =3K4 
0 ~0 


Now the mean value of cos? @ over the unit sphere is 4. Hence we have, approxi- 
mately, 


Kee | Tcos? 0 dw ~~ cos®6 | Ldwcilo 2 ee (2.47) 


This approximation is exact if J could be expanded in odd powers of cos@. Itisa 
well-known approximation (Eddington 1930) in the theory of radiative equilibrium. 
Hence (2.43) becomes 


aK /dr=—F[X(1—-e)+1). tit... (2.48) 
From (2.42), (2.47) and (2.48) we have a single differential equation for J, namely 
ay 
7at = ~ =3¢[X(1 —e)+1][J —40T4] 
or 
G4 Jade? = 8’c(.J to K*), | eee (2.49) 
if 
pen/{se[ Xb el ls ee een (2.50) 
The solution is 
J (7) =40T4+7-"(Ae-* Br+Be*Br), ss. eas (2.51) 


where A and B are constants of integration. 
Since J must remain finite near 7=0 we must have 


AtBales Gate: = bo (2.52) 
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To find A, we observe that, at the boundary [(«R,0)=0 for 6>7/2. Hence 
7/2 : 74] 2, 
F(«R)=27 i I(«R, 9) sin 6 cos 6 dé, J (aR) = 2n | I(aR, 6) sin 6 dé. 
0 


Replacing cos @ in the first integral by its mean value from 0 to 7/2, we find 


HERS OTGR)S «erreurs ere (253) 
But from (2.48) 
< aJ 
F(«R) = — alee (2.54) 


Hence from (2.51), (2.53) and (2.54), 
2aRT4 


[ we -§- 5] ore” )-[ peg §— 3 Jexr(—208) 


iS pnd for the surface flux F(«R)=oT*E,(aR) where es ae 
| B+2¢(1 ~ aR) ~ (2-2 E + aR | ) exp (—2pa2) | 
wet De . (2.56) 


Since the surface flux would be o7" if the boundary were replaced by a black 
_ body, E,(«R) represents the effective emissivity. It is clearly a generalization of 
the formula (2.11) for E(«R) to the case when « <1. 
Consider the value of E,(«R) for e=1 and «R large, i.e. condition E satisfied. 
We find 
Pico) =4e(B42eys. ®t es (2.57) 
Since 8 = 1/(3e) it follows that E,(co)=1-072; the correct value is, of course, 1. 
This shows that our approximation is quite close, and suggests that we write, for 
the effective emissivity, «,’ 
ental) 0-93 (an) ) ge) eS on ees (2.58) 
instead of the quantity L,(aR) defined by equation (2.56). 
The following special cases of formula (2.56) are of interest. 
1. For e=1 we should have, from formulae (2.58) and (2.57) «;'(aR) = E(«R). 
In Figure 2, curves 1 and 2 show how these quantities compare. Curve 1 is, of 
course, exact, and the two curves differ because of the approximations (2.47) and 
(2.53) which have been introduced to obtain expression (2.56). It will be seen 
that curve 2 lies just below 1 for most of the range. If £,(aR) were plotted it would 
almost lie on curve 1 except near the origin. ‘The curves show that the approxi- 
mation (2.58) is only a few per cent out. 
2. Suppose E is satisfied, that is «Ris large. ‘Then for uniform and Rayleigh 
scattering 
€no0) = 3°24 /el(2a/e+4/3). ewes (2.59) 
For diffuse reflection 
3:72/e (2.60) 


Eo ( 00) = Der a/i(13 —4e)3\. ocean. 


Figure 3 shows how these expressions depend upon e. 
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When E is satisfied, higher degrees of approximation are possible and ¢,'( 00) is 
better expressed by 
2(1—«) In {1+ /(3e)} (3e)F] 
———_ | 1- —— — |- -eee- 2.61 
5 | I aes ine 


This approximation gives e’=1 when « =1; it is almost indistinguishable from — 
(2.59) in Figure 3. This expression is derived in Appendix IT. 
3. When «R is small, that is, the cloud is in state “E’”’, E,(aR) =4eaR/3 and 


Hegre N gol = eee (2.62) 


As we should expect physically, the radiation is simply the total radiation 
emitted by each of the particles. 

The dependence of <’ upon a for a cloud of given size and mass is seen in Figure 
2, which is plotted for a cloud of aluminium oxide particles. The curves show 
several points of interest. There is not very much difference between the curves 
for diffuse and uniform scattering—the latter lying slightly above the former. 
For different values of «, the curves are of the same general shape, but reduced in 
scale along both the y and x directions. From the shape it will be seen that the 
radiation is sensitive to a, the particle size, over a fairly restricted range. Fore=1_ 
this range is approximately 10-3 to6x10-3cm. Whene=0-1 and 0-01 the corre- 
sponding ranges are 5 x 10-4 to 2 x 10-3 cm. and 10-3 to 3 x 10-*cm. respectively. 


€, (00) =1— 


2.3. Radiation from a Cloud of Particles in the Form of a Spherical Shell 


Suppose now that JN is constant within the spherical shell of bounding radii R’ 
and R but that there is no matter within the sphere of radius R’. The whole of the 
analysis of the radiation within the uniform sphere given above applies to the 
radiation within the shell, except the boundary condition (2.52). This is replaced 
by a boundary condition at y= R’. Since there is no matter inside the sphere of 
radius k’, there can be no net flow of radiation outwards across the surface. Hence, 
by symmetry, 

Fg) = 000 eae (2.63) 


‘The formulae for the effective emissivity corresponding to equation (2.56) is 
«(a R,R)=186 1-p{p+2«[ 1-2] 
: ‘ a 
1 aBR 1 —1 
+[exp{—20p(R-R)}]| He || B-2«(1+ =a) |} 
= p{p—2¢ [1+ | 
ap 


+[exp (2aP(R—R’)}] oo | | B +2 (1 = sel 


When k’ =0, it is easy to verify that (2.64) gives formula (2.58). 

Suppose a(—R’) is large. Then a fortiori «R is large and we find (2. 64) 
gives (2.59) or (2.60) for.the effective emissivity. For «R small we again obtain 
formula (2.62). If « is fixed, as R’— R, «,'(a, R, R')>0. But suppose we fix a 
and the total mass of Ae in the cloud but vary R’. As R’ increases, the 
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material is compressed more and more into a thin shell between radii R’ and R. 
| In fact, if « is. the value of « when R=0 the value of « in the general case is 


a oR (ReaRO Ie eee en ae (2.65) 


Figure 1. 
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Figure 2. Variation of effective emissivity with particle radius. 
Wik 3° /omnnk— Atta lne=— ll exacttormnuulal(2.13))502,,e— 1, 
approximate formula (2.58); 3, «=1, uniform scattering, 
uniform distribution in thin spherical shell, formula (2.66);. 
4, «=0-1, uniform scattering, formula (2.58); 5, «=0-l, 
diffuse scattering (2.58); 6, e=0-1, uniform scattering, 


Figure 5. thin shell (2.66) ; 7, e=0-01, uniform scattering (2.58). 
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2, uniform scattering, first approxi- 
mation (2.59); 3, diffuse scattering, 
first approximation (2.60); 4, uniform 
scattering, second approximation (2.61). 


In the limiting case R’>R we find «BR=aBR’— oo; exp{—2aB(R—R’')} 
exp {—28R/3}. Hence the effective emissivity of a thin spherical shell is 
given by 


3-72e(1 — exp { —2o8R/3}) 


colt RR) 59.4 (B= 2e)exp {— 2uqGRDY 


fe (266) 
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In Figure 2 we have plotted this expression against the particle size for e = 1 and 0-1 
and regular scattering. The emissivity is not markedly different from that of the 
uniform spherical distribution. This rather surprising indifference of «’ to the 
distribution law is also shown by calculating «,'(«, R, R’) as a function of R’. We 
have calculated e’ for particles of 0-004 cm. and a cloud of the same size and mass as 
for Figure2. The variation of <’ as R’ increases from zero to 1 ft. is only about 1%. 

The approximations (2.47) and (2.53) may not be very good as R’ approaches 
R but the general behaviour of (2.64) suggests that the radiation from the cloud is 
not affected very much by the density distribution of the particles within it. 


$3) COOLING-OF THE CLOUD: 
EXAMINATION OF ASSUMPTION (1) 

We must now examine the conditions under which assumption (1) is valid; 
that is, for the temperature to remain uniform as the cloud cools. Suppose that 
the particles are maintained in some way at a uniform temperature T for time 
t<t’. Att=v?' let the heat supply be discontinued so that the cloud begins to cool. 
Let G(r) denote the initial rate of fall in temperature. 


Gr== eich. |.) et ae (3.1) 


‘The condition that the temperature remains uniform in the subsequent cooling is, 
clearly, that G(r) is independent of r. 

Heat may be transferred from the inner to the outer layers of the cloud by 
radiation, conduction or convection. The condition that G(r) is independent of r 
will, in each of the three cases, impose restrictions on the physical quantities 
involved. We first consider the transfer of heat by radiation. 

The total loss of heat from a spherical part of the cloud of radius r equals the 
amount of heat emitted from the portion radius 7, less the heat absorbed from the 
outer shell of thickness (R—7r). Assuming that the density of the cloud is uniform 
this gives 

4nC |" G(r)Pdr=4aroTtE(ar)— A, es (3.2) 
where A denotes the heat absorbed from part (R—r) and C the thermal capacity of 
the cloud per unit volume. ‘To estimate A we consider as before the heat transfer 
in a thin cylinder SP’Q’T, in directions within solid angle dw about s. The por- 
tion P’Q’ lies in the sphere radius 7, and SP’, Q’T lie inthe shell (R—r). We 
have from (2.6) ; 
I?’,8)=(o Tn) Sexp(—aSP ie See (3:3) 


A fraction 1 — exp (—«P’Q’) of this light is absorbed in P’Q’. Hence the total 
amount of light absorbed by the inner sphere from the outer shell is 


A=8noT* |" [1 —exp(—aSP’)][1 —exp(—aP’Q)]z dz, 


where 2 is the distance of P’Q’ from the axis. ‘The integral cannot be evaluated, 
but the maximum value of SP is (R?—7?)? and the minimum value is (R—7r). 
Hence it follows that the heat absorbed lies between the quantities 

8roT*[1 —exp{—a(R—r)}] i a[1 —exp { —2a(r? — 2?)}}] dz 
and 

80 T4[1 —exp {—a(R?—r?)4}] |” s[1 —exp {—2a(r?— 2)}}] dz 


¢ 
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or 
| AnoT*?E(ar)[1 —exp {—a(R—-r)}] <A <4a0T*E(ar)[1 —exp { —a(R? —7?)#}]. 
aS (355) 
Bs ahetituting the value of A given by aa (3.2) we find 
r°E(ar) exp { — a(R? —1?)} <r aes G(r)r? dr <r°E(ar) exp {—a(R—r)}. 
See a (3.6) 


If G(r) =G, where G is independent of 7, we must have, for all values of r, 
rE(ar) exp { —a(R* —71?)8} <CG/30T* <r E(ar) exp { —a(R—r)}. 
The two outer limits r-+0 and r->R must converge, hence the condition that 
G(r) is independent of r is 
Lt rE(ar)exp{—a(R—r)}= Lt rE(ar) exp {-—a(R—7r)}. 
r—>0 . r—>R 


re (3.8) 
From (2.13) Lt r4£(ar)=4a/3. Hence (3.8) becomes 
r—>0 
do Rise exp (oI) (RR ee Se ey 6) i reret 8 (3:9) 
Now from equation (2.13) for E(x) we have 
Mees arp (a BNEG R= Dar (eye, | task (3.10) 
: 2 


where a, =[(n+3)(n+1)]/(n+2)!, m odd, and a,=n/(n+1)!, n even. Hence, 
since all the coefficients a,, are greater than 0, condition (3.10) is only satisfied if 


Cae Mleesvng ate Vr (FT) Oe TAPP leh S Ss 621) 
that is, if condition E’ holds. This is the reverse of E, the condition found for 
the radiation to be equivalent to that of a solid back-body radiator. 


Unfortunately relation (3.6) cannot tell us much about the behaviour of G(r) 
since it does not follow that 


ql kor) exp { — a(R? —7?)?}] — 7 G(r) or  [2E(ar) exp {—a(R—r)}]. 


In view of equation (3.6), however, and the fact that the above relation is certainly 
true near ay Z seems reasonable to suppose that the quantity 


G*(r)= = 507 dry ae TAOS {—a(R—r)}+exp {—a(R?—7?)*}], 0... i) 


gives a good approximation to the behaviour of G(r). 

In Figure 4, G*(r) is plotted as a function of r for the two cases «R =0-3 and 
aR=3. It willbe seen that in the first case G* is very nearly uniform but this is by 
no means correct for «R=3. It is easy to visualize the physical mechanism in 
cases aRS1. Strictly each part of the material radiates the same amount of 
energy but most of the energy is re-absorbed except near the surface where it can 
escape into outer space. Hence the temperature of the outer part is lower than 
that of the interior. ‘The same general considerations will obviously apply to 
particles for which e« is less than 1. 

We must now examine the possibility of the temperature being maintained 
uniform in state E by conduction of heat. Suppose heat is lost from the surface 
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of the cloud by radiation only, whereas inside it is transferred by conduction only. 
Let k be the thermal conductivity of the gas occupying the space between the 
particles. If H is the energy radiated by the cloud in unit time 


H=—4nkR*O0Tjdr\|_.0 > 7) |) eee (3.14) 
“R a 
But H==4n| cL ar. wer (3.15) 
0 6st ‘ 
If the temperature is nearly uniform, we may replace the second integral by — 
OT (® 3 47R3 OT 
An ib dra" Cae (3.16) 
where 07/0t represents the mean value of dT/dt. But we also have 
oT ar ror OF a or ; 
Dees aye hh Ry geal, 
ree | Catia or RC eee (3.17) 
Integrating, we find 
T(R) — T(0) =(CR?/6k)(OT/0t)=H/8rkR. —......... (3.18) 


For R=1ft. and T=2700°K. and a temperature difference of 100° c. between | 


the centre and surface of the cloud we should require k™1C.c.s. units. 


At T~2700°x., k for air is about 4 x 10~*c.c.s. unit; thus it appears that the — 


thermal conductivity is much too low to maintain uniform temperature. 
There remains transport of heat by convection or turbulence; this is difficult 


> 


to deal with quantitatively. Butit seems possible that this factor is quite important 


in view of the high temperatures involved. 


Hence we conclude that E’ is a necessary condition for the temperature of the — 


cloud to remain uniform if it were so initially, unless the temperature is maintained 
uniform by convection. 


§4. DISCUSSION 

Summarizing our results, we find: 

(i) If the transfer of heat by convection may be ignored, uniform temperature 
during cooling requires condition E’. When this condition is not satisfied, the 
cloud is cooler near the surface than in the interior. 

(ii) When the temperature is uniform and «=1 ,the cloud is equivalent to a 
solid black body if E holds. Whether or not E is satisfied for ¢ less than 1 the total 
radiation lies between the value it would have if « equalled 1, and this value 
multiplied by «. In other words, if e’ is the effective emissivity «<eé’ <1. 

Now let us consider the bearing of these results on Caldin’stheory. If transfer 
of heat by convection is important it follows that this method is legitimate provided 
that the appropriate expression for the effective emissivity is employed. On the 
other hand, if transfer of heat takes place through radiation only, it would seem at 
first sight that our results show that it would be impossible to apply Caldin’s 
theory to transient light sources, for this theory assumes uniform temperatures— 
which requires condition E’—and yet assumes E and that e =<’ in the determination 
of the heat radiated. But one or the other of E and E’ may be satisfied, not both 
simultaneously. ‘Thisisa little too drastic however ; for if E’ is satisfied, the theory 
is still valid provided that for the effective emissivity e’ of the radiator we use not 
the ordinary solid surface emissivity ¢, but the expression 


e-=3MeAak*dp,. eee (4.1) 
This follows from equation (2.16). 
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Condition E’ imposes a lower limit on the particle size—for example, with a 
‘4 cloud of radius one foot containing 4 grammes of magnesium oxide we require 
j a>2:0x10-3cm. Magnesium metal burning in air gives small cubic crystals of 
»| the oxide with side lengths from 10-> to 10-4cm., so that in this case E’ is apparently 
| not satisfied. But experimental evidence must decide the state of the cloud in 
9} each particular case. 

It is of interest, however, to consider some of the consequences if, both in the 
formation of the cloud and its subsequent cooling, the cloud is cooler near the 
| surface than in the interior. A necessary condition for this state of affairs would be 
the non-fulfilment of condition E’.. The amount of radiation given off is deter- 
mined mainly by the temperature of the surface layers. For, owing to absorption, 
1 the fraction of the emitted radiation which escapes to outer space decreases as we 
# move from surface to centre; defining the effective temperature 7, of the cloud by 
1 the equation 


HO) Oe 2 eee eg hn Nay (4.2) 

id then 7,<T. 

| In certain cases it is possible to find the effective temperature in terms of the 
9 temperatures in the interior and at the surface. For, assume E and that «=1; 
then equation (2.3) gives the intensity. Suppose the temperature, instead of 
_ being constant is given by 


Pee ecxpeails,  he (4.3) 


| This would represent a temperature distribution like that of Figure 4, curve 2, in 

1 general shape—nearly uniform temperature in the interior surrounded by a cooler 

layer. 7) is a rough measure of the thickness of this layer. The surface temper- 
MA 0 & ME Pp 


ature will be 
fb dey Bhd So 1 BO ees eis a (4.4) 


| if 7, is the temperature of the interior. The equation for J is 


cos 6 alr, 6)/07 =1(7, 6) —(¢T54/7)[1 — exp ( —7/79)]. 


ene (4.5) 
Integrating, for 0<7/2, 
LOD y= (67 7/7) bl —pro/(Zo+COSO)] 2.2 asics (4.6) 
Hence 
71/2 
F(0) =oT$=2n| 1(0, 6) sin 0 cos 6 dé, 
0 
=Glo(l —Duirgl —74ln(1-+75)/7o}]. 0 eee es (4.7) 


For 7,>1, that is for large optical depth of the cool layer, 7,=T,, and for 
| t)<1, or for small optical depth of the cool layer, T,=7). From Figure 4, 
curve 2, we should expect in practice ty~1, i.e. from (4.7) 


Pe SOT EEO OL TS a te (4.8) 


This result still holds even if the cool layer is so thin, compared with the radius of 
the cloud, that the mean temperature is effectively 7. 

Conditions E and E’ havea bearing upon the efficiency of transient light sources. 
Consider two clouds initially at the same uniform temperature, and of the same 
thermal capacity, one of which, S, satisfies E, the other, S’, E’. Let us also 
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assume that in S convection is not sufficient to maintain uniform temperature. 
After a long time both will have ultimately lost all their heat to the surroundings. 
But whereas S’ cools with uniform temperature, S has a cool outer layer and so its 
radiation escapes more slowly, and at a lower effective temperature than S. In 
the temperature range which we are concerned, the ratio of the quantity of light 
energy to total energy in the black body spectrum increases as the temperature 
increases. Hence, although S and S’ ultimately lose the same amount of energy, 
a greater fraction escapes as light in S’thanin S. Thus the efficiency of acloud asa _ 
source of light increases as we pass from state E to state E’. 

Another consequence of the cool outer layer is that the emergent radiation will 
not be black-body radiation for states E. When E is satisfied the intensity of 
radiation of wavelength A emerging at @ is 


2hc* (° exp(— sec ®@) ‘ 
1,0, 8) == | aise (4.9) 

This will not necessarily correspond to the Planck distribution at temperature 
T.. since T in the integrand is a function of 7. The energy distribution of (4.9) for 
radiative equilibrium has been studied by Milne (1900). In this case the spectrum 
of the total flux corresponds to a temperature 4% higher than T,. The spectral 
distribution is slightly bluer than black-body radiation for the effective temperature — 
T.. The deviation from the Planck formula is so slight that it is of little import- 
ance in astrophysics. It seems reasonable to conclude that the same applies to 
light sources. 


$5) ECAC. Se LHe ORY. OF ACOOL ING 


In view of the limitation on Caldin’s method it is natural to enquire if a more 
exact theory is possible in which account is taken of variations of temperature in 
the light source. The above analysis suggests various possible approximations ; 
for example we might use equation (4.8) to give the radiation lost in unit time 
assuming T,=0. It is quite easy, however, to set up an exact formal theory of 
cooling but the equation must eventually be solved numerically. 

Consider how the fundamental equation (2.27) is modified, if, instead of being 
maintained at constant temperature the material cools by radiation. The only 
change in the equation is that 7, instead of being constant, is a function of r. We 
have, also, an additional equation determining the rate of change of temperature at 
any point. 

For a single particle the total radiation absorbed from a beam of intensity J in 
directions within solid angle dw about directions s in unit time is ewa*Jdw. Inte- 
grating over all directions this becomes eza?J. ‘The amount of energy radiated if 
T is the temperature is 477a"«oT*, hence the equation to determine T is 


Spyh | 
aoe a =ena*J —4na%eoT* or a sel eT) hate ac (GU) 


where h is the specific heat of the solid material. 

If we assume a uniform density distribution, (2.27) becomes (2.35). The 
temperature and radiation intensity are determined for all ¢ and r from equations 
(5.1) and (2.35) provided the initial and boundary conditions areknown. (5.1) and 
(2.35) cannot be solved exactly. ‘The only feasible method seems to be to retain 
(2.35) unmodified, but to replace (5.1) by a difference equation. — 
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Suppose we know TJ as a function of r att¢=0. Then, solving (2.35) we can 


“4 determine I and so J at t=0. Suppose At is a small time interval, so small that J 
and 7 are not changed very much at the end of it. Then if AT denotes the small 
4 change in T in intervals At, we have from (5.1) 


NE AAP SAG TAZAT Aap oltre idiom sadek: (5.2) 


It only remains to express the solution of (2.35) in a suitable form, when T is 


a given as an arbitrary function of 7. Using the method of §2, (2.35) becomes 


Gr idm = Br —4oT sy aes (5.3) 
) since equation (2.49) is still true if Tis an arbitrary function of 7. 
Writing Dig espienp (Be) ledrs ~~~ ag (5.4) 
:: (5.3) becomes (dL/dr)+2BL+4BoT*rexp(—fr)=0. a... (5:5) 
Integration gives 
) L(x) =L(0)exp(—2fr) —4Bo exp (—28r) |" xT4(x) exp (Bx) dx, «0... (5.6) 
@ and for J 


J (7) =A exp (fr) 
+Bexp(-—fr)- ae exp (£7) ie exp (— 262) dz |" x xT*(x) exp (Bx) d 


‘ A and B are determined by the boundary condition at r=0 and r=R, namely 
F=0when7=0, F=2J whenz=za?NR. The total radiation from the sphere is 


Heer eye ae mt peak cida Phe. hose Grek (5.8) 


“There is no point in making any calculations by the above method until more is 


1 known about the size of the particles. 


AEP EAN Dihex. -T 
(a) Law of scattering for particles with a perfectly reflecting surface 
Let the angle the radius vector to any point on the surface makes with s be 0. 


i ‘The particle is placed in the path of a parallel beam and we require the amount of 
| light scattered into dw about s when the surface scatters light according to the law of 


reflection. Clearly all the light which meets the surface between radii whose angles 
6 lie between ¢/2 and ¢/2 + d¢/2 is scattered into directions s between angles ¢ and 


 d+dd¢. 
Hence, using the definition (2.25) of y(P,s’,s), we have 
(1/4:)y(P, s’, s)27 sin 6 dd = (sin? (f+ dd) /2—sin? g/2). ...... (6.1) 
Hence VERRSe Ste be me ME Socices (6 2) 


and thé scattering is uniform. 


(b) Law of scattering for particles with a diffusely reflecting surface 
Let OT (Figure 5) be normal to the plane containing s ands’. OR lies in the 


plane containing s, s’ and OQ lies in the plane containing OT and OR. 
: 9-2 
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All the points T, Q and R lie on the surface. Put 7 ROS=0, 7SOS=0, 
ZROQ=p. Let yy, #y be the angles between the normal at Q and the directions 
sands’. Now the element of area at Q is a?coswdudé. 

If J is the intensity of the incident light in directions in dw about s, the amount 
of light falling on this area at Q is [cos ,dwa*?cosydud6. By Lambert’s law of 
diffuse reflection, the amount of light scattered into directions dw’ about s’ is. 


a 1dw' cos pI cos bs, dwa? cos pd dd. 


But cos#,=cospcos6, cos%,=cospucos(?—8s). Hence the above expression 
becomes 77! dw dw’a*I cos* 4 cos 6 cos (8 ~8)dudé. Hence 


1 oe? Ure 
—y(P,s’,s -=] { os? cos 8 cos (6 —5) dé d, 
rest! Vii) aan RRA eR LO 


= alla 8) coso+sid|.  «- «. _-— ee (6.3) 


If we put ¢ equal to the angel between directions s and s, i.e. 6 =z —8, then we have 
finally from (6.3) ; 
y(P,s’,s)=(8/37)(sind—dcosd). sew (6.4) 


APPENDIX II 
Calculation of the higher approximation to « for the case when E is satisfied 


When « =1 we have from (2.8) 
I(S,s)=(cT4/7)[l1—exp(—aPS)]  —_........ (7.1) 

. if P is the point where direction —s through S cuts the surface of the envelope of 
the cloud. If condition E holds, «PS is large except when PS<R. Hence J 
is uniform and of magnitude o7“/m except in a thin layer near the surface of the 
sphere. At the surface J cannot be uniform; for when 6>7/2, J=0. As the 
thickness of this layer is reduced, conditions near the surface clearly approximate 
more and more to those of a plane surface. Hence with E satisfied, the second 
term of (2.33) is small compared with the first and 


QI/as=cosO(AI/ar). nae (7.2) 


The same kind of behaviour is to be expected for the case «<1. Hence (2.27) 
becomes 


al(r, 8) eo l™ 1 ai, ; ; 
csp =natN | — I(r, 0) + San I(r, 8’)y(r, 6, 8) de | 
oe (7.3) 
or if Tre [,Nar, VE oA (7.4) 
Ol(7, 8) 1—e)/ , ; : fh 
cos 6) = I(r, 0) — oa I(z, 8’)y(7, 8", 0) dex’ — — Paste ee: (7.5) 


When 7 =0, r= R and from (7.4), 7 is large for points in the interior of the cloud 
not near the surface. 

An exact solution of (7.5) is possible only for the case y(r, 6’, 0) =1, that is, for 
scattering by specular reflection. Equation (7.5) can then be transformed into an 
integral equation which is very similar to the fundamental integral equation of 
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the theory of radiative equilibrium developed by Schwarzschild and others for 
astrophysical purposes (Milne 1928, Titchmarsh 1937). The solution is compli- 
cated, however, and requires extensive computations (Lindblad 1923). 

The method of §2 can, however, be applied (with only slight modifications) to 
obtain a first approximation for <’. From this approximation it is possible to 
calculate a second approximation in the case y=1. The method cannot be 
applied to the more general equation (2.35). For J, the first approximation is 


IAG st Bexp(— Ur), © 86 Ss (7.6) 
where f is given by (2.50) and 
Boe Lada LCA ) el a Se neeti i ths Libs ayes (7.7) 


Substituting expression (7.6) for J in (7.5) the equation for J is 


Cla we Le (1- 22 -j 
COSA at Vartipae aah er expC Bz), ieee (7.8) 


Integrating with respect to 7 from 0 to «, 
oT Bile) 


1. 0) — * 4n(1-+ B cos) if 0 Tr) 2, ee eees (7.9) 
Hence 
a : 2(1—e) In(1+ 8) 
= ee a ew: EE |S [poe ee ca 
F(0)=2n | 1(0, 8) sin @ cos 0d0 =oT E ao: {1 : i 
CE (7.10) 
It follows that if we denote the second approximation to e’ by [e,/(7)] 1 
; = 2(1—e) In {1 + +/(3e)} 
[ee (©) =1—- Jet af GO) 1— a, Pogue (7.11) 


For «=1 we find [e,'(0)];,;=1. 
In Figure 2 the first and second approximations are plotted. It will be seen 
that there is very close agreement between the two, suggesting that (2.59) and 


(2.60) are very near to the correct values. For diffuse scattering the method 


fails unless further approximations are introduced. It is interesting to note that a 
more sophisticated method of obtaining successive approximations to the solution 
of the fundamental equation of radiative equilibrium, which has been developed 
recently by Chandrasekhar (1944), could also be applied to equation (7.5). Since, 
however, it is not our aim to obtain very exact estimates of « we shall not pursue the 
matter further. 
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On the Schumann-Runge O, Bands emitted at Atmospheric 
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ABSTRACT. The Schumann-Runge bands of O,, previously analysed in the region 
3100-4200 a. and photographed below 3100 a. under small dispersion, with only a tentative 
analysis, have been obtained from 2000 a. to 5000 a. on sufficient dispersion for the vibra- 
tional anaiysis to be extended and for the rotational structure to be clearly resolved- 
Conditions for their emission are discussed and the probable process, involving dissociation 
followed by recombination, is outlined. On this basis the results obtained in discharges, 
ozonizers and high-voltage arcs are explained. Flory’s suggestion of predissociation in the 
upper state is reviewed and the experimental evidence found to be against it. ‘The conclu- 
sions reached are applied to considerations regarding the spectrum of the night sky. 


Silo IINTADIROIDNGHE API Kd) 

ONSIDERABLE interest is now being taken by spectroscopists in the spectra 
of the aurorae and the night sky, since from these a great deal of information 
may be obtained about processes occurring in the upper atmosphere. In 

connection with this work it is necessary to study the spectra emitted by the 
atmospheric gases in the laboratory. The present paper describes some new 
experiments on the emission spectrum of Og, and it is convenient to begin with a 
brief summary of previous investigations. 


(a) In the red and infra-red regions of the spectrum absorption bands corre- 
sponding to the transition a‘L;<-x*°X, have been observed. Near these bands 
others are found resulting from the transition 'A,<x*X,. Recently Kaplan 
(1947) has obtained some of these bands in emission from an oxygen afterglow. 

(b) Absorption bands found by Hopfield and others in the far ultraviolet 
region have been shown to be due to transitions from the O, ground state to 
various Oj states, forming a Rydberg series of bands. 

(c) Hopfield (1930) discovered a number of emission bands near 2000. 
forming a single v” progression. ‘They are produced in a condensed discharge 
in a O,-He mixture and are attributed to O, on experimental evidence only. 

(d) Herzberg’s ultraviolet absorption system (°,<-x°X, ) is forbidden and 
has not yet been obtained in emission in the laboratory, though some strong night- 
sky bands are thought to belong to this system. 

(e) The most readily observed absorption in O, is that found by Schumann, 
which sets the limit to observations in the ultraviolet region with air-filled spectro- 
graphs, ‘These bands result from the transition ?X;<-x*X, and bands observed 
by Runge in a high voltage arc in the region 4200-3100 a. are known to be part of 
the same system in emission. 

(f) At very high pressures new absorption bands appear in the oxygen 
spectrum, and these are thought to be due to the molecule O,. Further absorption 
bands are obtained in solid and liquid oxygen, most of which have been analysed 
into four systems. 
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The present paper gives extended data on Runge’s bands and discusses the 
processes involved in their emission. The new results are not consistent with the 
suggested predissociation of the 3X; state at v' =2, and, therefore, all the relevant 
experimental evidence is examined in order to settle this question definitely. 
Finally the conclusions reached in this investigation are reviewed in relation to 
the work on the night-sky spectrum. 


§2, EXPERIMENTAL 


The discovery of the Schumann-Runge bands in emission was published by 
Runge (1921), and various workers measured and analysed the plates he obtained, 
the most complete analysis being given by Lochte-Holtgreven and Dieke (1929). 
This work was confined to the region 3100. to 4200 a., and no spectra were pub- 
lished. In his experiments Runge used a high voltage D.c. arc in flowing oxygen. 
Whilst investigating the oxonizer discharge, Wulf and Melvin (1939) compared its 
spectrum when filled with cylinder oxygen with that given by a high voltage, 
low current a.c. arc in the same oxygen. They published a small dispersion 
spectrogram of the two sources; that of the arc showing Schumann-Runge bands 
extended from 4000 a. to 2500 a., whereas that of the ozonizer shows only N, 
Second Positive bands. Millon and Herman (1944) have carried out experiments 
with discharges in extremely pure oxygen and also report the appearance of some 
of the bands; this work will be considered in a latter section. 

In the present investigation a high voltage, low current a.c. arc was used. A 
two-litre Pyrex flask contained the source, which consisted of an arc running 
between two thin (0-03 in. diameter) platinum-rhodium wires. The Pt-Rh wires 
were supported on the ends of tungsten wires sealed into ground-glass joints which 
fitted into the flask, so that easy adjustment of the electrodes was possible. A 
Hyvac pump and a McLeod gauge were provided and the gas was supplied from a 
cylinder. Pressures between 0-05 mm. and 1 atmosphere were employed. ‘The 
discharge is a true arc, the electrodes becoming red- or white-hot and presumably 
maintaining the arc by thermionic emission. Various electrodes were tried, the 
difficulty being to find a material which can remain red-hot in oxygen without 
oxidizing. Platinum has been found to work reasonably well, but it is an advantage 
to use a platinum—rhodium alloy as the melting point is then raised. 

The voltage in such a source is determined mainly by the separation of the 
electrodes, whilst the current can be varied by means of variable resistances in the 
primary circuit of the transformer used. Usually the arc ran at about 1000v. 
and 0:-lamp. These figures are, however, only of qualitative significance since 
the voltage characteristic is far from sinusoidal. Current and voltage wave-forms 
were investigated with a Cossor double-beam oscilloscope and showed the form 
predicted for a.c. arc (Simon 1905, Loeb 1939, figure 286); the current 1s practi- 
cally sinusoidal, whereas the voltage shows a sharp breakdown peak at the start of 
each half-cycle, followed by a period of constant voltage. Higher intensities 
could presumably be obtained by using higher currents, butthe source is convenient 
in its present form, runnning extremely steadily ard giving sufficient intensity for 
photographs to be obtained on a Quartz Littrow (Hilger E.1) spectrograph with 
exposures of about } hour in the more intense regions. 

Commercial cylinder oxygen was used and, although this contains various 
impurities, especially nitrogen, no bands of molecules other than O, appeared on 
the plates taken at 1 atmosphere pressure. The bands are in fact quite strongly 


The Schumann-Runge O, Bands PLZ 


Set Pee MESS 1O NS PROCESS 

The present experiments show that the intensity of the band system falls 
rapidly with decrease of pressure in the discharge, so that at } atmosphere it is 
exceedingly weak. This behaviour is completely different from that of many other 
gases (for example N.) which emit best in a discharge at low pressure. Hence it 
becomes necessary to attempt an explanation of this result. The reason for this 
behaviour may best be discussed in conjunction with the figure, which shows the 
potential curves of the upper and lower states of the Schumann-Runge system. 
These curves are drawn to agree with the experimental evidence, absorption, 
emission and photochemical. As is well known (e.g. Mott and Massey 1933), 
a collision of an electron with a molecule leading to electronic excitation may be 
treated according to the Franck-Condon principle. Thus it follows from the 


° 


oA) +00) 


Potential (Volts) 


7 (Angstrom Units) 


‘figure that the excitation of O, molecules in the ground state will most frequently 
result in dissociation in the upper state. The experiments of Ladenburg and 
Van Voorhis (1933) on the absorption coefficient of oxygen in the vacuum ultra- 
violet region show that the transition of maximum probability (equal instantaneous 
‘internuclear distance in the upper and lower stat s) gives an average kinetic energy 
of about 1-7 ev. per atom after dissociation. The action of a discharge on oxygen 
will thus be to start dissociating it, and various authors have pointed out that 
atomic oxygen is very reactive, readily combining with gaseous or solid impurities. 
However, in the high pressure experiments here described there is a much greater 
probabuity than in the low pressure discharge that the atoms will recombine again 
to the upper state of the Schumann-Runge system. ‘This would probably take 
place by a three-body process, and the molecule would then return to the ground 
‘state by emitting the Schumann-Runge bands. It should be noted that the dis- 
-sociation leads to a normal O(?P) atom and an excited O(1D) atom. But since the 
1D level is metastable (mean life 100 sec. (Condon 1934)), it will not easily revert 
to the ground state by radiation, and so prevent the recombination taking place. 
Furthermore, recent experiments (Vegard and Kvifte 1947) have shown that the 
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de-activation of metastable oxygen atoms by collisions with other oxygen atoms. 
or molecules has a very low probability. The high rotational temperature of the 
bands is consistent with these views of the emission process, since the recombin- 
ation of atoms with high kinetic energy will obviously lead to the formation of 
molecules of high rotational energy. 

It is of interest to apply the above reasoning to the results obtained by Wulf 
and Melvin (1939) with an ozonizer discharge in commercial cylinder oxygen. 
At 1 atmosphere pressure they obtained no oxygen bands but only N, impurity 
bands. Here the dissociation phenomenon gives an explanation, since, although 
in the ozonizer the volume of gas in which excitation takes place is much larger than 
in the case of the high voltage arc, the amount of energy dissipated may not be very 
different (Wulf and Melvin used the same transformer for both). Hence to a first. 
approximation the intensity of bands excited by direct cdllisions in each case will 
be the same. But in the case of oxygen recombination followed by emission is 
favoured in the arc with its small volume of highly dissociated gas rather than in the 
uniformly distributed dissociation of the ozonizer. 

Millon and Herman (1944) have run discharges in extremely pure oxygen 
during investigations on the forbidden red and green lines of oxygen. At a 
pressure of a few centimetres Hg the Schumann-Runge bands occur with intensity 
comparable to that of the red and green lines. It is difficult to tell from their short 
paper the absolute intensity of the bands, but since the lines simultaneously 
observed are forbidden it seems safe to assume a low intensity. ‘This then falls 
into line with other results since in any but the purest oxygen impurity bands. 
would obscure the O, bands and the oxygen atoms would tend to react with the 
impurities rather than reform O,. Millon and Herman’s wavelengths of heads are 
not of sufficient accuracy for a completely certain vibrational analysis of the bands, 
but apart from the measurements of Lochte-Holtgreven and Dieke, the assign- 
ments (0,11), (1,11), (0,10), (1,10), (0,9) probably agree with the present work. 
However, the bands they identify as (0,8), (0,7), (0,6), (0,5) are possibly the ones 
taken here as (2,9), (2,8), (2,7), (2,6), which seems more consistent with the expected. 
Condon parabola. They have also found a band at 2331 a. which they take as. 
(0,4); no band at this wavelength has been observed in the present investigation. 
It is hoped that a complete rotational analysis will settle any doubts about the 
vibrational scheme. ; 

When the major part of this paper had been written a note by Lalji Lal (1948). 
appeared. He has observed bands from 2450 a. to 4490 a. in a high frequency 
discharge in pure oxygen. ‘I'wenty-six heads have been measured on small. 
dispersion, and these bands have been assigned to the Schumann-Runge system. 
The analysis he gives is quite different from that given here, and no bands at his. 
wavelengths can be found on the present spectra. ‘The two results have been 
briefly compared by the author (Feast 1948), and further data on Lalji Lal’s bands, 
especially a photograph, are needed before the matter is completely discussed.. 


§). PREDESSOCTLAT LON 


Flory (1936) has suggested that the behaviour of oxygen can only be success- 
fully explained by assuming that the upper state of the Schumann-Runge system is. 
predissociated by a repulsive °I1,, state at v' =2 of the ®X) state, this leading to the 
formation of two normal oxygen atoms. This would mean that bands with v’ >2: 
would not be observed in emission. However, according to the analysis proposed. 
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here, bands with v’ =3 are found, and it seems difficult to find any other interpre-- 
tation for the unanalysed region at the short wavelength end of the spectrum than 
that it is due to bands with v'>2. It is therefore important to re-examine Flory’s 
arguments to see if all oxygen phenomena can be explained without assuming: 
predissociation. 


(i) Emission Bands 


Flory suggests that from Lochte-Holtgreven and Dieke’s results a Condon 
parabola may be deduced which indicates that the bands (3,15), (3,16), (3,20), 
(3,21), (3,22) should have appeared on Runge’s plates: in fact these bands are not 
found. He maintains that collision quenching could not have prevented obser- 
vation of the bands since he assumes Runge’s pressure as less than 0-01 atmosphere.. 
This can be criticized on two grounds: (a) As has been shown earlier, the bands are 
only produced with high intensity at high pressures ; Runge’s pressure would thus 
probably be near atmospheric, where the collision frequency is about 101°sec.,. 
while the radiative lifetime is about 10-$sec. Hence collision broadening is. 
very important at these pressures and explains the tendency to broadening in the 
short wavelength bands. (b) The intensity of the bands, on the plates published: 
here, indicates a wider Condon parabola than that taken by Flory, so that the bands. 
quoted above do not lie on the line of maximum intensities. 


(ii) Absorption Bands 


Analysis of the absorption Schumann-Runge bands has been carried out by 
Curry and Herzberg (1934) working at 1 atmosphere pressure and by Knauss and 
Ballard (1935) at 5-15mm. Hg. Curry and Herzberg’s photograph is thought by" 
Flory to show broadening for v’ >2, thus bearing out the predissociation theory. 
The effect must in any case have been small to escape the experimenters, and the 
fact that in the lower pressure work the bands are clearly resolved, with no trace of | 
diffuseness, suggests that any broadening may be put down to collision quenching. 


(iii) Fluorescence Experiments 


Rasetti (1929) found O, fluorescence at atmospheric pressure but none at 
8mm. Hg. This is explainable on a theory of predissociation, but it may also be 
understood if collision quenching does not fall off more rapidly than the number of 
molecules excited by fluorescence when the pressure is reduced. 


(iv) Photochemical Experiments 


The pressure decrease when oxygen is exposed to radiation of A<1751 A. is. 
due to the reaction of the atoms formed by absorption of the Schumann-Runge 
continuum with wall impurities (Neujmin and Popov 1934). Flory and Johnston 
(1935) used 1 cm. of air as a filter and claim that the small pressure decrease, still 
noted, is due to a predissociation following absorption of radiation of A>1751 a. 
However, since the observed effect is small, there are some other factors worth 
considering: (a) 1cm. of air may still allow transmission of a small amount of 
radiation with A<1751a., thus forming atoms in the normal way; (4) though, as 
Flory states, the effect of dissociation by a bimolecular collision of an excited and. 
a normal molecule will be small, it may not be completely negligible ; (c) absorp-- 
tion of radiation in the region of the Herzberg system continuum is a possible: 
mechanism for dissociation at wavelengths greater than 1751 a. 
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(v) Conclusion 


From the above discussion, and from the results on the emission bands, it 
appears that the experimental evidence is against the predissociation theory. In 
certain cases the results previously quoted in support of the predissociation are 
found to have other possible causes. 


$6. APPLICATION TO THE NIGHI-SKY —SP ECR EM 


Since oxygen is abundant in the upper atmosphere, its presence or absence in 

‘the spectrum of the night sky is of great importance. At present there seems a 
-good deal of agreement amongst astrophysicists that the Herzberg oxygen bands 
occur in emission. The presence of the Schumann-Runge bands is somewhat 
doubtful, though various authors claim to have identified them. The Herzberg 
emission bands are calculated from the absorption measurements and from the 
-ground-state differences of the Schumann-Runge bands (Av”). At present the 
Av" values required have to be extrapolated from the measurements of Lochte- 
Holtgreven and Dieke. The work described here was started with the object of 
extending the data on the Schumann-Runge bands so that the Herzberg bands 
could be calculated without the uncertainty of an extrapolation. Since it appears 
that the heads are too weak for direct measurement, the necessary information will 
-only be obtained by a complete rotational analysis, and it is hoped to publish this 
shortly. Various workers (e.g. Pearse 1943) have calculated the wavelengths of 
bands corresponding to bands of the Schumann-Runge system with v’ = —1 and 
—2 and have noted some good coincidences with night-sky bands. ‘This has led 
to some speculation as to whether the accepted wv’ values are correct. None of 
these bands can be found on the present spectra, and hence it must be assumed that 
the v’ assignment is correct. The most conclusive test would be a study of the 
isotope shifts. 

We know that above 100 km. the oxygen is almost completely dissociated by 
solar radiation during the day (Wulf and Deming 1938). ‘Thus it seems possible 
that a three-body recombination to the upper state of the Schumann-Runge bands 
will take place at night, though this would yield low intensities. 

Swings (1947) and Nicolet (1948) have advocated that night-sky bands should 
be identified by first computing the low temperature profiles of the bands. In this 
connection it should be noted that if the Schumann-Runge bands do occur due to 
the process outlined above then their rotational temperature may be considerably 
higher than the accepted night-sky value, and the profile computing method would 
break down. 

Finally, since the pressure is low in the upper atmosphere, collision quenching 
will not be important and higher v’ values should be expected. 
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ABSTRACT. Wavelengths and intensities of twenty-one CN tail bands, observed in 
the carbon arc spectrum, are listed and v’, v’’ values assigned. Nine of the bands have been 
previously reported as occurring in the nitrogen afterglow excitation. The intensities of 
the tail bands are shown to be in accordance with the accepted potential curves for this. 


system. 


§1. INTRODUCTION 
AN is well known, the CN violet band system (22 —7%) occurs readily in 


many laboratory and astrophysical light sources. The system consists of 

several strong sequences degraded to the shorter wavelengths, together 
with some weaker bands degraded in the opposite direction and known as the tail 
bands. The principal sources used by investigators of this system have been the 
carbon arc in air and the nitrogen afterglow with carbon compounds present. 
In active nitrogen both the main bands (Jevons 1926) and the tail bands (Jenkins 
1928) have been studied. Extensive rotational structure measurements of the 
main bands as obtained in the arc have also been made (see Jevons (1926) for 
references). However, there seem to be no data available on the arc excitation of 
the tail bands. Such data are of importance both for laboratory identification of 
the bands, the carbon arc being used extensively in analysis, and for identification 
of astrophysical spectra. The present communication gives wavelengths and 
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intensities of twenty-one tail bands observed in the arc spectrum: nine of these are 
bands which have been reported in the afterglow spectrum, whilst the rest corre- 
-spond to previously unreported transitions. The band intensities are in qualita- 
tive agreement with the Condon parabola deduced from the accepted potential 
curves. 
§ 2. MEASUREMENTS AND ANALYSIS 

The measurements of the band heads were made from Hilger E.1 spectrograms 
-of a 110v. 5amp. carbon (graphite) arc in air. To obtain the tail bands with 
sufficient intensity it was necessary to have the main CN bands considerably 
-over-exposed. | 

The new data are given in the table. 


Ueno" Ac (A.) I An (A-) ve(cem™) (cm) %e—¥% 
Ss AS a 4078-72 

14, 14 4028-38 e 6? 4029-33 24817 24779 m 38 
{Lei 3 8984-52,6) 9 3984-62 25090 25049 m 41 
Wile 3944:52e 9 3944-66 25345 25288 m a7 
fils Ul 3909-90 e 10 25569 25498 m 71 


b 

10, 10 a b 

ial a b 
1 10 305)/:007em eS 3658-13 27332 DID fran 57 
LOS HOMES GU b Uf 30) 27477 m 75 
Ores 3602:63 e 7 b 27750 27642 m 108 

ee lel a b 
1210 3465733 6 b 28849 28818 m 31 
i 3432292 6 3433-00 29126 29080 m 46 
tO ears 3404-84 7 29362 29302 60 
ONE, 3380-29 6 29575 29494 81 
Soa 3359-06 4 29762 29651 111 
Tes) 3340-56 4 29926 29776 150 
6, 4 6322-29 4 30091 29868 223 
Dae: 3296-29 3 30328 29932 396 
10; 7 3203-46 4 31207 31154 53 
YE 3180-22 4 31436 Splat 65 
oh SH De 4 31637 31556 81 
7, 4 Sa 257 4 31812 31708 104 
O53 3127-56 3 31964 31827 137 
he 7 SEAS 2 32101 S929 182 


v’, v'’, vibrational quantum numbers; ., J., wave length, intensity,of band head in arc; 
Ay, wavelength of band head in afterglow spectrum (see text); ve, wave number correspond- 
ang to Ae; vv, the wave number of origin (‘“‘m”’ denotes measured values); ‘‘a’’, band 
-obscured by main bands; “‘b”’, band, but not head, observed in afterglow; ‘“‘e’’, measure- 
ment of band made difficult by strong overlapping main bands. 


In calculating vg, the combined data of Jevons and Jenkins have been used; 
where the experimental values are available they are indicated by ‘‘m’’. For the 
majority of the bands reported here the values of B’ and B” are not known with 
‘sufficient accuracy for the positions of the band heads to be reliably calculated and 
compared with the experimental values. However, the values of v, — vp given in the 
‘table are all of the correct order and vary from band to band in the expected manner. 

The column A, in the table gives the wavelengths of heads measured by Jenkins 
in the afterglow spectrum and it will be noted that in several cases there is a 
considerable difference between this value and X,._ This is partly due to the fact 
‘that in the arc many of the tail bands are badly confused with the main CN bands 
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(those most affected are marked “e”’ in the table). A second reason for this 
discrepancy is that A, is the wavelength of the last rotational line in the R branch, 
as measured by Jenkins, the rotational temperature being such that no lines on the 
returning part of the branch are observed. A consideration of Jenkins’ measure- 
ments shows that the line taken as forming the band head may be 0:5 cm“! or more 
to the long wavelength side of the true head. 


§3. DISTRIBUTION OF VIBRATIONAL INTENSITIES 

The figure shows the distribution of the intensities in the CN main and tail 
bands as excited in the arc. The band intensities are head estimates, and since 
V¢— Vo 1s not constant from band to band, these values are not reliable for the over- 
all intensity of the band though they may be used for a qualitative investigation. 
The main band intensities were obtained by microphotometry of Hilger E.1 
spectrograms (‘Tawde 1936). No attempt has been made to put the two parts of 
the system on a-common intensity scale as the difference of intensities between the 
main and tail bands is too great to be estimated by eye. ‘The tail band intensities 
are eye estimates. 


yr 


The Condon parabola in the figure is drawn from the potential curves calculated 
by Herzberg (1939, figure 177) and is practically the same as that obtained from 
the Morse and Rydberg potential curves given by T'awde. It is seen that it fits 
quite closely to the experimental line of maxima. Bands near the parabola, in the 
(1, 0) sequence from v’ =4 to 8 and in the (0, 0) sequence from v’ =5 to 10, are not 
observed since these bands are either headless or fall in regions where the main 
bands are very strong. 

The vibrational intensity distribution in the afterglow (Jenkins 1928) is similar 
to that given here except that bands with the head far from the origin appear 
relatively weaker, due to the low intensity of high J value lines. 


§4. NOTE ON THE (4,4) MAIN BAND 


The wavelength of the (4,4) band is 3850-95..; this measurement does not 
appear to have been recorded previously. 


124 S. P. Sinha 


ACKNOWLEDGMENTS 


The author wishes to express his thanks to Dr. R. W. B. Pearse, who suggested 
this investigation, for his interest and advice. The work was made possible by a 
D.S.I.R. grant. 

REFERENCES 
HERzBERG, G., 1939, Molecular Spectra and Molecular Structure, Vol. 1 (New York : 
Prentice-Hall). 
JENKINS, F. A., 1928, Phys. Rev., 31, 539. 
Jevons, W., 1926, Proc. Roy. Soc. A, 112, 407. 
Tawpe, N. R., 1936, Proc. Indian Acad. Sct., 3A, 140. 


Ultra-Violet Bands of Na, 


ByiS. Py SINHA 


Science College, Patna, India * 
MS. received 27th Fuly 1948 


ABSTRACT. The ultra-violet bands of Na, between 43600 a. and 4 3200 a. have been: 
photographed in absorption in the first order of a 21 ft. concave grating with a dispersion. 
of about 1-3 a/mm. and their vibrational analysis carried out. All the bands belong to a 
single system and can be represented by the equation 

v=29342+119-33 (u’)—0:53 (w’)?—159-23 (u’’)+0-726 (u’’)?+-0-0027 (u’’)3, 
where u=v+}. 

The heat of dissociation for the upper state is calculated to be 6510 cm~?. The bands: 
are due to transition from the ground !33 state to an excited state which is probably 137 
in character. The dissociation products for the upper state are either 3 ?S+4?P or 
3 *S+3 ?D atoms of sodium. 


§1. INTRODUCTION 


been studied by Walter and Barratt (1928), Weizel and Kulp (1930), 

Kimura and Uchida (1932) and Sinha (1947). In an earlier report (Sinha 
1947), it has been considered that the bands belong to seven different systems of 
which only three appear at all extensive, the remaining ones being merely frag- 
mentary. ‘Two of these systems, reported there as system 1 and system 3, have 
been found to be strong while system 2 is comparatively weak. The analysis. 
proposed suffers from the defect that the vibrational differences are not quite 
satisfactory. It was, therefore, considered necessary to photograph the bands at 
higher dispersion and determine the molecular constants more accurately. This. 
has, however, been done only for the bands lying between A3500 and 3200a., 
corresponding to the first ultra-violet system mentioned above. 


Te ultra-violet bands of Nag lie in the region 43650-A 2500 a. and have 


§2. EXPERIMENTAL 


The bands have been studied in absorption. The details of the absorption 
tube used have been described elsewhere (Sinha 1948). Sodium, freed from the 
oil in which it was stored, was placed at the centre of the tube and was heated 
electrically from outside. ‘To prevent rapid distillation of the metal to cooler 


* The work described in the present paper was done at Imperial College, London. 
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parts some nitrogen was introduced inside the chamber. The bands were found 
to appear most satisfactorily when the temperature inside the tube was about 
750° c. and the pressure of nitrogen was 5 to 10cm. of mercury. A hydrogen 
discharge tube of the type described by Hunter and Pearse (1936) was used as a 
source of continuum and ‘the bands were photographed in the first order of a 
21-ft. concave grating with a dispersion of about 1-3 a. per mm. 


§3. APPEARANCE OF THE SPECTRUM AND MEASUREMENT 


The Na, ultra-violet bands first appear at about 700° c. when the pressure of 
nitrogen inside the chamber is about 5 cm. of mercury and lie between A 3450 and 
A3250 a. As the temperature is raised they become stronger and extend up to 
A 3600 a. on the long wavelength side and 43200 a. on the short wavelength side. 
With temperature higher than 850° c., the absorption in this region becomes 
almost continuous and quite strong bands begin to appear at wavelengths lower 
than about A3000 a. At still higher temperature (1000 to 1100° c. or more), the 
entire region from A 3600 to A2800 a. is completely absorbed except for a narrow 
region at about 43100. which still continues to be transparent. 

The A 3600—A 3200 a. system of bands which forms the subject of the present _ 
investigation appears most satisfactorily when the temperature is about 750° c. 
The strong bands lie up to A3550 a. while fainter ones extend to even longer than 
43600 a. Closely spaced lines pertaining to rotational fine structure appear over 
the entire region, but as the structure of one band overlaps that of the neighbouring 
band in most cases, it is not possible to separate the structure of each individual 
band, nor is it possible to see the band origins. A rough calculation with values 
of B’s of the two states involved reveals that the head is formed by as low as the 
third or fourth member of the head forming R branch, and since both B’ and B’ 
are so small that the low member lines of each branch lie very close together, the 
individual lines near the band origin cannot be distinctly resolved, with the 
consequence that the band origin cannot be observed. Even the band head, 
although quite distinct, was not as prominent as is ordinarily expected in view of 
the close spacing of the lines near the head. ‘The reason for this is that in the 
rotational structure of a band there is an intensity maximum which lies at aJ value 
that is higher the higher the temperature or the smaller B. In the present case 
B is low and the temperature used was high, and therefore the higher members of 
the branch were stronger than the lower members which formed the heads. ‘The 
difficulty could probably be overcome by working at a lower temperature and a 
correspondingly longer absorption path. 

Measurements are given in Table 1. The iron arc was used for comparison 
spectra. The values of )’s given here differ from those published before (Sinha 
1947). Since the heads are formed from relatively weak lines there is difficulty 
in correctly locating the position of the band-heads in both cases, but since the 
present measurements are from a plate of much higher dispersion and since the 
bands were more satisfactorily developed it is considered that the present measure- 
ments are definitely the more accurate. 

The bands were all degraded to the red and appeared to have single heads 
suggesting 1Z-1Z transition. It may, however, be noted that even if the bands 
had more than one head, it may not always be possible, under the circumstances 
described above, to see the different heads distinctly and the band may appear to 
have just a single head. 
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Table 1. Na, Bands between A3600 and A3200 4. 
d / ur Yobs— ir , 7a 
lO) Int. v,U Aait (4.) Int. v,v 
Veale 

3595-2 1 Daw 2 3400-4 8 ghas 
86-7 1 — — 3395-9 5 PS | 
80-2 1 ewe 2 92-0 4 1,0 
77-0 1 21 3 87-1 5 Ae 
70:5 1 0, 9 3 82-4 6 Seed 
65-4 2 4, 12 g 77°8 4 F.0 
62:3 3 oud. ) 74-1 5 ae) 
59-0 2 2, 10 2 69-2 10 4, 1 
55-5 3 1-9 4 64-6 4 50) 
52-0 2 Ca 2 61-4 5 6, 2 
Appiah Ws 4,11 2 56:5 8 set 
44:3 3 3, 10 2 51:5 a. 4, 0 
40-8 3 eo 4 48-7 5 9 ee 
37°3 2 1c.8 2 43-7 6 Gs vi 
33-5 2 O. 7 3 38-8 10 5) 
30-1 2 4,10 1 31-4 5 7a 
26-3 3 Cus 3 26-3 10 6, 0 
22-6 3 aes D 19-2 5 Oe 
19-0 3 REG 2 14-0 10 7, 0 
15-0 3 O06 3 07-3 5 9, 1 
12-1 3 479 3 3295-5 5 10, 1 
08-3 3 ses 1 90-0 9 9, 0 
04-5 3 Ce) 3 83-9 5 aoe 
00-6 3 t 6 2 78-4 8 10, 0 
3496-6 3 Os 2 72:3 4 tad 
94-2 3 7a) 3 66:8 7 ican) 
90-2 3 ny) 3 66-2 5 14.0.9 
86°3 4 Da 2 60-8 4 135k 
82:2 3 eae 3 12.70 
78-2 3 0, 4 2 23 ? 15, 2 
72:2 3 BONG 2 50-2 6 Ie 
68-1 5 eax’ 2 49-6 3 14, 1 
64-0 3 hy 28 1 44.5 3 16, 2 
60-0 3 0, 3 0 39-8 5 foe 
50-0 7 3. 4 2 38-7 3 jieug 
46-0 4 pars 0 34-1 3 17Lre2 
41:6 4 ue 0 33-2 5 14, 0 
De) 7 046% 0 29-5 5 1oaa3 
27-7 5 eD 0 28-2 3 164 
23-2 4 0, 4 0 23-6 3 (Py 
18°5 4 303 isi 19:3 4 2003 
14-0 5 Dee 0 17-6 4 17,04 
09-4 4 ie 0 13-4 3 OLE 
; 0, 0 0 09-2 2 pate 3 
Oe : aa 4 | 03-3 2 20, 2 


Yobs— 


Veale 


Ste 2 ose Orso Shiga S ofee Moor Mee ocr ater) Glee 4 ery 0.9$t 0. gr ee 
£.0r1 L.IPY £.yhr 9: SHI o.LF1 0-6F1 €.0S1 O.1S1 9.1S1 0.SS1 


Z.QS1 g-L4S1 ‘UIp uvsyay a 
Ae ISIle RAL 
bSOLE cH 6OZIE OZ 
Ss 0.86 96 ot 36 
< 9S60E HH TLILE 6L 
S.96 66 A 96 
LS80€ ch ELOTE SI 
6.66 66 D oor 
SSLOE % CLONE KH 69OTE LI 
S.1or EON Sion 
LISO€ & 8960E OT 
Goaton ZOl Sron 
60L0E & L980E ST 
0.ZOI IOI oe Lor 
8090€ ch F9LOE GH OZ6OE +1 
0.gor gor oa 
a 8S90E el 
8 o.Lor Lot 
= ISSO€  6OLOE 76)! 
— S.Lo1 ; gol C2 it 
= EtbO€ & TO90E LT 
a S.Lor Lol © gor 
S DECOE ch FOFOE Ol 
in 0. gor gol oo gol 
RQ S7ZOE cA OSEOE 6 
nig 0.601 601 ce 
~~ 6110 x 8 
2 O.OII OIL 
NS bS86Z cH GOODE & 99TOE L 
1 Scale CAs SOT Ta Sr arate 
= ZHL6T HA 86867 ch SSOOE 9 
ss Sateiese on CN aii Sie Sar 
67967 % F8L6~ cH TH667 S 
~ o.€11 VII Ee ZIl 2 €11 
6087 & GLISZ\E OZEST E SOPS~ HX 1198Z ZIECT ch STS6Z % ZLIGZ M 67867 $ 
g-Sr1 Oi © OL Sie VED GN ested eee yaler QI, oo Grr git gil 
€COLE & F90D8T F OTST H OSEST K 96F8Z EK EF98T FZ ZLEZ brlOT ch OOHOHZT H 9SS6Z A ETL6T e 
Z.QII OD Si GiT ae NOR SVORT SO grr A ort 9 Se git Or peOwCi En wenl i: pot 
LORAZ  SPOLT H OOOSC F FETC E OSEST R L7SBZ EH 9L9BT A YET A LLOST H STLOT HH E~OT HR GEHOT EH 16562 @ 
€.L11 e LI CSAS) ST Ru | XS fais Oo fri aa ER dares, Ea seh Nia hie beer 
LII8Z £ OTST HR GOPEZ E SSS8Z cH GOLET A 9ST GH LIO6T H 99T6T A ZZE6Z 2 O8r6Z I 
Boles CLT eC TNS ST Oke Ta Dorks eee SET eee git "gil i QI gII 
. 666L7 & Sh18Z H COC8C § THH8C | 16S8Z Ht CHLBT A EOSBT BHT H 40767 A TIL6T 0 
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§4. VIBRATIONAL ANALYSIS AND MOLECULAR CONSTANTS 


The vibrational quantum numbers assigned to the bands are given in column 3 
-of Table 1. Table 2 gives their wave-numbers in vacuum in the usual wv’, v” 
square array. The asterisk denotes the position of the Na, band that was over- 
lapped by the A3303 a. line of Na. The mean differences between the vibrational 
levels for the ground state from the work of Loomis and Nusbaum (1932) are given 
at the foot of Table 2. 

Nearly all the bands measured can be represented by the equation 

vy = 29342 + 119-33(w’) —0-53(u’)? — 159-23(w”) +0-726(u")? + 0-0027(u")?, 

where u=v+4. The terms in uw” correspond to those given by Loomis (1928). 
The values for v,,.-Veare given in column 4 of Table 1 have been calculated from 
the aboveequation. Although there is not much discrepancy between the observed 
and the calculated values there seems to be rather a systematic difference for large 
values of both the quantum numbers v’ and v”. The discrepancy for large values 
of v" seems to be due to fact that the constants in v” used in the above formula do 
not refer to the band heads but to the magnetic rotation lines which very nearly 
represent the band origins, whereas the measurements included in Table 1 refer 
to the band heads. The difference between v,,,,4 and Voigin for bands degraded to 
the red becomes larger for higher values of the vibrational quantum numbers of 
the lower state (Jevons 1932), and since for such bands v.44 — Vorigin 18 Positive, 
Vops — Veale Should also be positive. ‘This agrees with the present observation. 
The discrepancy for large values of v’ suggests that higher terms in v’ should be 
included in the above equation. It further indicates, as is often seen with the 
excited states of such molecules, that the levels are tending to converge more 
rapidly than is demanded by the linear relationship between w, and v. However, 
the discrepancies are nowhere too large and the arrangement seems satisfactory. 

Vibrational constants for the four states of Na, are given in Table 3. The 
heat of dissociation for the upper state of the present system has been calculated 
by the method of extrapolation due to Birge and Sponer and would probably be a 
little higher than the true value. The constants v,, w,, @,.%,) @,Ve are given in 
cm! and are those in the usual expansion 1(v)=v,+w.u—w.xu?+w,yu, 
where u=v+H4. 


Table 3. Vibrational Constants for the Ground and the Excited States of Na, 


Dissociation 


State Ye @e Wee @eVe » Products 
De 0 159-23 0-726 —(0-0027 6150 BRS a aS 
1 14680-4 TATE6 0:38 -- 8100 ayssy SpE 
ly 20318-9 123-79 0:63 —0-0094 2840 Betsy Suele 
Wr (2) 29342 119-33 0:53 _ 6510 Sosy 83 
or 
Seas) abep 


§5. DISSOCIATION PRODUCTS 


The energy possessed by the dissociated atoms in the upper state of the system 
can be calculated from the equation v,,5=¥%,9 +D’—D", where the terms used 
have the usual significance. Since, v» 9=29360cm"', D’=6510cm"! and 
D’ =6150 cm", we have v.44, =29720cm~!. Also from the line spectra data it is 


—_ 
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known that for sodium 32S—32D=29170cm™ and 32S —4?2P =30270cm™!. 
The above are the nearest atomic states to the calculated v,,,.,, other states being 
still farther apart. Since, however, the energy possessed by the dissociated atoms 
is midway between that corresponding to 32S —32D and 32S —42P, it is difficult 
to decide between these two sets of dissociation products. Of the three terms 
used in calculating v,,,,,,D”, the heat of dissociation for the ground state, is known 
correct up to +175 cm™ due to the work of Loomis and Nusbaum (1932). Table 2 
giving the v’ —v” scheme indicates that the value of vp 9 can be taken to be quite 
reliable. Hence the only uncertainty lies in the value of D’ employed in the above 
calculation. D' has been calculated by Birge and Sponer’s extrapolation method. 
This method is applicable to an electronic state for which the relationship between 
@, and v is linear and has been found to yield fairly satisfactory results for the 
ground states of non-polar molecules. At any rate, the-smaller the range of 
extrapolation employed the more certain one can be of the calculated value of the 
heat of dissociation. For the excited states the values calculated by this method 
are often too high. A critical discussion on the suitability of this method has been 
provided by Gaydon (1946). Although in the case of the present system of bands, 
over the range investigated the relationship between w,, and v seems fairly linear, 
it is difficult to predict how it would behave for higher values of v’. Judging from 
experience of other states of such molecules it can be concluded that, possibly, the 
true D’ and hence »,,,,, may be smaller than the value given above. This would 
lead to the conclusion that the dissociation products are 32S +3?D atoms rather 
than 37S +4?P atoms, but a definite decision should not be based on this con- 
sideration alone. 
§6. DISCUSSION 


The ground state of Na, is a stable 1X; state formed from two atoms of sodium 
each in the unexcited 37S state. The other state formed from these two atoms 
is a°X state which is regarded asarepulsive one. Since in absorption, transition 
is to take place from the ground 1%f state, the only excited states that matter are 
those of the type 1X and HII. The excited states are obtained by considering 
one of the sodium atoms to be raised from the ground 375 state to higher and 
higher excited states like the 32P, 42S, 37D, 4?P etc. If we consider states 
resulting from 3 2S and 3 ?P atoms, there will be as many as eight of them of which 
only one will be 12 and only one 4JI,. ‘There should thus be two systems of 
bands 1X <-1Z; and 4I1,<'2; whose upper state should dissociate into 37S 
and 32P atoms of sodium. These correspond respectively to the red system 
studied notably by Fredrickson and Watson (1927) and Fredrickson and Stannard 
(1933), and to the green system studied notably by Loomis (1931) and Loomis and 
Nusbaum (1932). 

Let us now consider that one of the sodium atoms is raised to the next excited 
state i.e. the 42S state. The 32S+4?S atoms of sodium will give four states, 
out of which there will be one 1X, and no 1II, state. ‘The transition from the 
ground state to this 12;* state should give rise to a system of bands whose origin 
should lie at about A 3900 a., if we assume that the excited state has nearly the same 
dissociation energy as the ground state, or at longer or shorter wavelengths 
according as D’ is greater or smaller than D". Experiment does not reveal any 
band due to Na, near 43900a. The bands at longer wavelengths are the green 
system of Na, bands which are definitely known to be due to 3°S+3°P atoms. 
The bands at shorter wavelengths are the bands under present investigation and 
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lie between 43600 and A3200 a. If we wish to associate this system of bands with 
Na 328+42S, D' for the upper state should be about 2500cm™?. Although 
this is not an improbable value, the results of analysis do not suggest that D’ 
could be as low as that. Further from analogy with the known & states of Nap, 
K, and Nak, it is expected that this state should also be fairly deep. ‘This might 
lead one to conclude that the 12*(32S+2?S)<1L; system of bands is not 
observed. 

The next excited states should result from consideration of an unexcited 375 
and an excited 32D atom of sodium. This gives twelve molecular states, with 
only two possible transitions, viz. 1X and 'IL, from the ground state.- Conse- 
quently two systems of bands should be obtained whose upper states will dissociate 
into 32S +3?D atoms. Again, just as two systems of bands are observed due to 
transitions from the ground state to the excited states resulting from 37S +3?P 
atoms of sodium, similarly two systems should be observed due to 1X; <-1L; and 
1]],>12X; transitions whose upper states will dissociate into 37S +4?P atoms. 
We should thus have four systems of bands, two associated with 32D atoms and 
two with 4?P atoms, all of which should lie close together. Actually, however, 
only three systems of bands are observed in the corresponding region including 
the present system (Sinha 1947). Study of the absorption spectrum of K, has 
also revealed only three systems of bands in the corresponding region (Sinha 1948). 
Work is in progress to explain this discrepancy and also to understand the absence 
of bands associated with 325 +4?5S atoms. 


ACKNOWLEDGMENTS 


The author wishes to acknowledge his indebtedness to Assistant Professor 
R. W. B. Pearse, D.Sc., for his kind help and valuable suggestions, and to Patna 
University, India, for the grant of their Birla scholarship. 


REFERENCES 


FREDRICKSON, W. R., and STANNARD, C. R., 1933, Phys. Rev., 44, 633. 

FREDRICKSON, W. R., and Watson, W. W., 1927, Phys. Rev., 30, 429. 

Gaybon, A. G., 1946, Proc. Phys. Soc., 58, 525. 

Hunter, A., and Pearse, R. W. B., 1936, ¥. Sci. Instrum., 13, 403. 

Jevons, W., 1932, Report on Band Spectra of Diatomic Molecules (London: Physical Society), 
p. 55: 

Kimura, M., and Ucuipa, Y., 1932, Sci. Pap. Inst. Phys. Chem. Res. Tokyo, 18, 109. 

Loomis, F. W., 1928, Phys. Rev., 31, 323. 

Loomis, F. W., and Nussavum, R. E., 1932, Phys. Rev., 40, 380. 

SINHA, S. P., 1947, Proc. Phys. Soc., 59, 610; 1948, Ibid., 60, 436. 

Wa ter, J. M., and Barratt, S., 1928, Proc. Roy. Soc. A, 119, 26s. 

WEIZEL, W., and Kup, M., 1930, Ann. Phys., Lpz., 4, 971. 


| 
| 


131 


On the Equations of Motion of Crystal Dislocations 


By F. C. FRANK 
H. H. Wills Physical Laboratory, University of Bristol 


MS. received 25th August 1948, and in amended form 25th November 1948 


ABSTRACT. It is shown that when a Burgers screw dislocation moves with velocity v 
it suffers a longitudinal contraction by the factor (1—v?/c?)}, where c is the velocity of 
transverse sound. ‘The total energy of the moving dislocation is given by the formula 
E=E,/(1—v?/c?)t, where Ep is the potential energy of the dislocation at rest. Taylor 
dislocations behave in a qualitatively similar manner, complicated by the fact that both 
longitudinal and transverse displacements and sound velocities are involved. 


§1 

RENKEL and Kontorowa (1938), studying the properties of what is essentially 
} a one-dimensional dislocation model, found equations of motion for the 
dislocation which can be brought into a form strikingly analogous to those of 
a particle in special relativity. It has a limiting velocity equal to the characteristic 
velocity of sound in the (actually dispersive) system, suffers a contraction analogous 
to that of Lorentz as it approaches that velocity, and the ratio of its energy in motion 
to its energy at rest likewise depends “‘relativistically’’ on velocity. This result 
(which we actually discovered independently) will be discussed more fully, along 
with other properties of one-dimensional dislocation models, in papers to be 
published with Mr. J. H. van der Merwe. It has been noted before (Frank 1948) 
that in the three-dimensional case of crystal dislocations the kinetic energy of a 
Burgers ‘screw’? dislocation (Burgers 1939), extrapolated by the approximate 
formula valid at low velocity up to the velocity of transverse sound waves, is equal 
to half its potential energy. This is equivalent to the relativistic equation 
my = Ey /c?, where Ey is rest energy and my rest mass. It is shown below that all the 
elementary relativistic formulae are applicable to a screw dislocation, and that the 
behaviour of a Taylor “‘bridge”’ (or ‘“‘edge’’) dislocation (‘Taylor 1934) may be 
described as “relativistic with complications” : the limiting velocity is that of 
longitudinal sound waves, but the velocity of transverse sound waves also plays a 

part. 

SF 

The displacement function u(x, y, 2) of aresting dislocation according to linear 
isotropic elastic theory, given by Burgers (1939) for bridge and screw dislocations, 
is a solution of the three equations (1,, 1,, 1,) of elasto-static equilibrium typified 


by 
Ge GeO 0? 0? GC" 4 
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oy 
Dislocations in motion must be represented by solutions of the three equations of 
elasto-dynamic equilibrium (wave equations) typified by 


0 (ou Ou, | Ou 0? 02 02 au, 
ea tua aa tag tas bir (VAT Sigh et eas 
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In particular, a dislocation in steady motion with velocity v in the x direction must 
be represented by a propagational solution of equations (2), a time-independent 
function of x’, y and z where x =(x—v#). With this substitution, x is replaced by 
x’ and (02/022) is replaced by v?(0?/0x’?) in equations (2), producing three equations 
(3,,, 3,, 3.) which are equivalent to equations (1) with x replaced by x’, except for 
the operands of (02/0x’?) which are : 


in (3,) : (A+2u—v%p)u, in place of (A+2p)u,, 
in (3,) : (u—wvp)u, in place of pu, 
in (3.): (u—v%p)u, in place of pu,. 


§3 
Now one of the solutions of (1) is the screw dislocation along the z-axis, for 


which (Burgers 1939) 
Ue, u, =0, U,=\b)/ Zar) ave Can WjiX)y mes (4) 


equations (1,) and (1,) reducing to 0=0 and (1,) to pd?u,/dx? + we?u./dy? =0. 
Likewise for a screw dislocation in motion, equations (3,) and (3,,) reduce to 
0 =0 and (3,) to (uw —v%p)07u,/0x'? + nd?u,/dy2=0. Hence in this case equations (3) 
are brought into identical form with equations (1) on introducing 


x" =2'/(1 — o%p|p)! =(x —vf)|(1 —0*/e2)}, 


where c, =(z/p)? is the velocity of transverse sound waves in the material. Hence 
the displacement function of a screw dislocation in motion is identical with that of a 


> 


screw dislocation at rest, apart from a “‘ Lorentz contraction’’. 


§4 
The potential energy of a resting screw dislocation is 
By = 4p § [Ou,/@x)? + (u,|y)?] dx dy dz = (ub2/47) In (14/70) 


per unit length in the z direction. 19, of the order of the interatomic spacing 8, 
represents the boundary within which linear elastic theory is inapplicable, and r, 
the outer boundary of the material. Moderate changes in the size or form of 
these boundaries are unimportant. Neglect of strain energy in the non-linear 
elastic region within 7 is justified when 7,/7, is large; but on the evidence of the 
one-dimensional model it is likely that the results obtained below will remain true 
in a higher approximation which takes account of this energy. The potential 
energy E, of a moving screw dislocation is given by the same integral, which on 

writing dx =(1—f?)! dx", where B=v/c, becomes 


E, = 344 | [(Ou,/Ax")"/(1 — B®) + (Au,/0y)*](1 — B*)* dx" dy dz 
= 2B o/(1 — B?)* + 3Eo(1 — *)* = Ey(1 — $8?)/(1 — B?)*, 


since u is the same function of («”, y, z) as it was formerly of (x, y, 2). 
The kinetic energy associated with a moving screw dislocation is 


Ey =4p | (0u,{0t)* dx dy dz =p | (u,[0x")?[02|(1 — B2)|(1 — B®) dx" dy de 
= {0*/(1 — B°)!}(p[u)Bo|2= 3 BoB2/(1 — BP) 


On the Equations of Motion of Crystal Dislocations 1a 


Hence the total energy associated with the moving dislocation is 
- B=E, + E,=Ey/(t— "Jt. 


For dislocations moving at low velocity we may define an equivalent mass m, 
by equating the energy of motion E—E, to }mpv?. Then m= F,/c?. 

$5 

The work done by an applied shear stress o on the glide plane when the 
dislocation is translated a distance x is (o . b)x per unit length of dislocation. This 
energy goes partly into kinetic energy and partly into setting up oscillations in the 
crystal. This latter “radiation damping” of the motion of the dislocation should 
arise from the acceleration of the dislocation, under the influence of the applied 
stress and in its interactions with irregularities in the lattice, e.g. other dislocations, 
and also, at high velocity, from the periodic structure of the lattice itself. Formally 
we may define the force acting on the dislocation, f, as (o.b) diminished by a 
frictional force arising from the radiation damping. We expect this correction to 
be small in a reasonably perfect crystal, provided that the velocities or accelerations 
are not exceedingly high. ; 

We could identify the energy of a moving dislocation by integrating energy in 
the material. (It should be noted, however, that E— Ej, which might reasonably 
be called the kinetic energy of the dislocation, is not quite identical with the kinetic 
energy in the material, whichis £,.) ‘The same cannot be done with momentum, 
but we can define an equivalent momentum, p= Jfdt, corresponding to 
E-—E,={fdx. Then introducing v/c = 8 =sin¢ (a substitution which simplifies 
all the formulae of elementary relativistic mechanics) so that E=£,/cos 4d, 
we have dp/dt =f =dE/dx =(Eysin ¢/cos? ¢) df/dx =(E,/c cos? f) dd/dt, and hence 
p= J (E,/c cos? d) dd =(Ey/c) tand =(E/c?)v and E? = EG + p’c? as usual. 

A screw dislocation (unlike a bridge dislocation) may move in any plane 
containing its singular line. The transformations employed above may be 
performed independently for motions in the x and y directions. Hence p has 
independent x and y components, i.e., isa vector. is ofcourse scalar. In view 
of the “radiation damping” we cannot suppose there is strict conservation of the 
components of p unless we are able to assign a corresponding equivalent momentum 
to oscillations in the crystal. 

$6 

The case of the Taylor bridge dislocation is more complicated, since both 
longitudinal and transverse displacements occur in its motion. If it lies along the 
g-axis, with displacement vector b on the x-axis, which must also be its direction 


of motion, the displacements u, and all derivatives with respect to z vanish by 
symmetry. Equations (1) for the dislocation at rest then reduce to 


(A +211) 0740,,/0x? + 20?u,,/Oy? + (A+ p)O?u,/oxdy=0, (5,) 
(A+ 2p)0?u,,/Oy? + p0?u,/0x% + (A+ )Ou,/Oxdy=0. ss... (55) 
The major displacements are w,, parallel to b, and to a first approximation 
satisfy (5,,) with the final term neglected. ‘Through the final term of (5, the 
x-displacements give rise to y-displacements which react back upon the 
x-displacements through the final term of (5,)- This description is appropriate 
to a solution of (5) by successive approximations. 
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The solution of the first approximation is 


U4) = (6/2) arctan [{(A+2p)u}4(y/x)]. wee (6,): 

The closed solution obtained by Burgers (1939) is 
u,,=(b/27)[O+ {((Atp)/(A+2)}sin@cos6], nee 7) 
ay = —(b/2m)L{ue[(A-+2p2)} In (r/7o) (A+ m)(A+2p)}sin29], 0... (7,) 


where tan@=y/x and r=(x? + y?). 

The difference between wu) and uw, as given by (7) is small, being at the worst 
0:014v (when 6=7/3) in the case A=p (Poisson’s ratio c=}). Thus the first 
approximation is a not unsatisfactory one. 

Equations (3) for the bridge dislocation in motion become 


(A+ 2u —0%p)02u,,/Ox"? + pd2w,,/Oy? + (A+ y)O2u,/Ox’dy=0, ..- (8,,) 
(A-+2)02u,/Ay? + (u — v%p)02u, (Ox? + (A+ u)O'u,/Ox'dy=0. ....-. (8,) 


To the same first approximation as we have just discussed, taking the first two: 
terms of (8,) the problem reduces ‘relativistically”’ to the static problem on 
substituting «” =x’ /(1—v?/c?)? where c?={(A+2y)/p}?, the speed of longitudinal 
sound waves. 

Note added in amendment 24th November 1948.—Further incomplete treatment 
of the problem (which originally followed here) can be omitted, as a complete 
solution has now been obtained, and will shortly be published, by Mr. J. D. 
Eshelby. From his treatment it appears that the requisite differential equations 
can be satisfied both by an irrotational and a distortional dislocation with displace- 
ment vectors in the direction of travel, or by any combination of the two. Both 
suffer ‘“‘ Lorentz contraction”’, with limiting velocities equal to the velocities of 
longitudinal and transverse sound waves respectively. Either motion alone can 
only be maintained if a force, normal to the direction of motion, is applied at the 
centre of the dislocation, with counteracting forces on the remote surfaces. These 
forces do no work, but there is no means of applying the force at the centre. This. 
force vanishes when the irrotational and distortional dislocations are combined in 
appropriate proportions (dependent on the velocity). This, therefore, gives the 
only physically important solution. ‘The force at the centre cannot be made to: 
vanish for a dislocation travelling faster than transverse sound. 
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LETTERS TO THE-EDITOR 


The Resistivity of Thin Metallic Films 


Various formulae have been suggested from time to time to account for the variation of 
the resistivity of thin metallic films with the thickness of the film. The simplest of these is. 
an empirical expression due to Planck (1914) : 


P=Ppo (1+ A/D) 


where p is the resistivity and D the thickness of the film, and p, the resistivity of the bulk 
metal; A is an unspecified constant. This formula accounted satisfactorily for his own 
experimental results and for Pogany’s classical work on the resistivities and optical constants 
of thin films (1916). 

When the thickness of the film is a small multiple or submultiple of the length of the 
mean free path of the electrons the number of electronic collisions is going to be increased.. 
The Figure illustrates the shortening of the mean free path / and the corresponding increase 


Direction of ae 
Electric Field ese 
—> ee 


in the collisions. AB represents the mean free path in the bulk metal, and, in this case, the: 
number of collisions will be v(1+Al/D) where k is a numerical factor. As a rough first 
approximation the latter may be said to be equal to 4/7 since 2/7 is the mean of the possible 
values of sin 8, and there are two surfaces. When D is larger than J, k will have a different 
(smaller) value, which can be calculated by suitable averaging. Substituting the above value 
for the frequency of collision between the atoms and electrons in Drude’s classical formula 
it follows that 


Ne? z On 
2mv(1--ki/D) °~ T442nD 
where a is the conductivity of the film, D its thickness, o~ the conductivity of the bulk 


metal, N the number of free electrons per cm*, and e and m the electronic mass and charge 
respectively. From (1) we obtain 


O==Poltt AlimD) 6 we at ptistews (2) 
which can be identified with Planck’s empirical formula mentioned above. 


Considering small temperature changes at about room temperature it is assumed that 
po and / experience appreciable, but D negligible, variations. Then (2) can be written 


p=opa(aT+a)(1+4lo/7DyT). Scee 205 (3) 


a is the normal temperature coefficient of resistivity, y the coefficient of the mean free path, 
a that part of the resistivity which does not depend on the temperature; the suffix (9) refers 
to values at 0° c. but T is measured on the absolute scale. 
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When (3) is differentiated with respect to T an expression is obtained for the effective 
temperature coefficient a’ : 


1 1 
(pee | Arecncors Gt 
pe (era) @) 
It is readily seen that a’ is zero when 
ja=T(1l-aDyT/4l), a ee etre (5) 


and for correspondingly smaller thicknesses or lower temperatures a’ will be negatives 
In general, a’ is less than a. 

It follows that when /, is large the thickness D for which a’=0 will be large for given 
values of yand T. There is evidence suggesting that the mean free path of the electrons in 
bismuth is exceptionally long. It is thus to be expected that comparatively thick bismuth 
films will exhibit a negative temperature coefficient (Tulley 1946). In addition, it follows 
from (5) that the thicker the film the lower the temperature at which a’=0 (Curtiss 1921). 
It is not suggested that expression (3) has a general applicability at lower temperatures, 
indeed, it is known that the resistivity varies at temperatures lower than 0° c. with a higher 
power of TJ'than unity. It is none the less interesting to note that the results of van Itterbeek 
and De Greve (1945, 1946) on the resistivity of nickel films at comparatively low temperatures 
are in good agreement with the above expressions. ‘The temperature coefficient of resistivity 
varies with the thickness, and the minima of resistivity vary with the thickness and the 
temperature as predicted. These authors quote the results of De Haas and Van den Berg 
(1937) on gold wires in support.of their own. Since gold wires cannot really be compared 
with thin films a somewhat different explanation is required. This will be dealt with in the 
near future. 


South-West Essex Technical College, R. A. WEALE. 
London E.17. 
November 1948. 
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Note on the Slowing Down of Mesons 


Fermi, Teller and Weisskopf (1947), and Fermi and Teller (1947), have given a formula 
for the rate of loss of energy.of slow mesons passing through a metal, the metal being 
treated as a degenerate electron gas. The formula is (Fermi and Teller, equation (15)) 

dW 2mm erV2 1370 

= = In ————————“———_ ln ee ee (1) 

dt 3a h c 
where m, e are the mass and charge of an electron, V the velocity of the meson, W its 
energy and vp the maximum velocity of the electrons which form the degenerate Fermi gas 
in the metal. ‘The formula is derived for the case when* V<vp. Some doubt has been 
thrown on the validity of this result (cf. for example Frohlich 1948); the purpose of this note 
is to confirm the result by showing that it is identical with that obtained in a problem on the 
electrical resistance of alloys, for which the application of the same methods has given 
results in agreement with experiment. 

If the effect of the lattice field is neglected, the calculations required to obtain the 
following two quantities are identical: (a) The rate dW/dt of loss of energy by a heavy 
charged particle moving with velocity V through the gas. (b) The rate of loss of energy of 
an electron gas moving with velocity V past a stationary charged particle (e.g. a proton). 
"The second problem is the same as that of finding the residual resistance of a metal containing 


* For a treatment of the effect of free electrons when V>vy see Kramers (1947). 
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i) dissolved hydrogen. For suppose that the resistance of a metal containing xN stationary 
| protons per unit volume is xp; Nis here the number of metal atoms per unit volume. Then. 
) ! the 1ate of loss of energy per unit volume, if a current of density j is flowing, is 


C0] eee a ns ey OP ee Ne (2) 
| If the metal is monovalent and thus contains N electrons per unit volume, we may set 
=~ j=NeV, so that (2) becomes xpN?e?V?. The rate of loss of energy due to a single stationary 
d| charged particle, which can be equated to dW/dt, is given by 


dVadi= nee Zettel oe =a yee | eatin. (3) 


The resistance p due to any type of foreign atom can be calculated by the method. 
originally due to Nordheim (1931) (cf. also Mott and Jones 1936, p. 289). We may write 


p=mv,A/e’, o aeskitons, / (4) 
where A is the area defined by 


AE a (OA cosa: # OG Wits eA 2,2 (5) 


I(8) dw is here the cross-section for scattering of an electron having velocity vg by a foreign: 
atom into the solid angle dw. This can be calculated if the field round the foreign atom is. 
known; in Mott and Jones (1936), p. 97, a screened Coulomb field of potential 


WG) (e2/7) exp. (9h aa eet TE, (6) 


y is taken, and q obtained by using the Thomas-Fermi method to find the density of the- 
electron gas round a proton. One finds 


gine ie ONi ata (ope: he Oi atiy. er (7) 


With this field [(@) can be obtained, using the Born approximation; one finds (using the- 
value given for 1(4) by Mott and Jones (1936), p. 294, equation (86))* from (4), (5), (6), (7) 


27e? thug 
O= In —). 
mMvy° e 


Substituting in (2), using the equation mu=(37)'rh/N*, we obtain the formula 
dW 2 mel? (=) 


ad 37 FR mN ege 


which differs from that of Fermi and Teller only in the 7 in the logarithm. Since the Born 
approximation is strictly valid only if Avg > e®, the change is unimportant. 

No experimental results are known for the increase of resistance due to hydrogen in 
monovalent metals, to which alone this analysis should be applicable; but the theory agrees 
well with experiment when applied to calculate the increase of resistance due to the addition 
of 1% of Ni, Zn, Ga, Ge or As to copper (Mott and Jones 1936, p. 293). The resistance is 
found experimentally to increase in the ratio of the numbers 1, 1, 4, 9, 16 for these metals. 
This can be explained (Mott and Jones 1936, p. 292) quantitatively on the above theory, 
taking a field round each foreign atom of potential (ze?/r) exp (—qr), where z takes the 
values —1, 1, 2,3 etc. It thus appears that, for example, a zinc atom replacing a copper 
atom behaves like an additional charge e, and the residual resistance is that due to a screened 
charge. 


H. H. Wills Physical Laboratory, ; N. F. Morr.. 
University of Bristol. 
23rd November, 1948. 
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REVIEWS OF BOOKS 


Preparation and Characteristics of Solid Luminescent Materials, by G. R. FONDA 
and F. Serrz (Editors). Pp. xv +459. (New York: John Wiley, 1948.) 30s. 


This book comprises twenty-nine papers together with associated discussions which 
formed the subject matter of a Symposium on Luminescence held at Cornell University 
in October 1946. It is concerned with recent researches on specific aspects of luminescence 
in synthesized phosphors, such as sulphides, silicates and tungstates. The papers presented 
at the Conference have been arranged in five main groups, while a further group contains the 
final discussions of the Conference. 

In the first group of papers a general survey of phosphors, their characteristics and 
methods of preparation is made. The next group contains four papers which deal with 
two important aspects of recent research on the storage of energy in phosphors. ‘These 
aspects are the mechanisms of phosphorescence and thermoluminescence, both of which 
arise from the thermal ejection of trapped electrons, and the processes involved in the 
infra-red stimulation of phosphors. The third group of papers is devoted to factors 
affecting the fluorescence characteristics of phosphors, such as phosphor constitution, 
nature and intensity of exciting radiation and the phosphor temperature. The next six 
papers deal in more detail with the storage of energy in phosphors and its release by thermal 
or optical ejection of trapped electrons. 

The last group of papers covers miscellaneous aspects of luminescence, such as photo- 
chemical action in phosphors, the function of manganese as a luminescence activator in 
different solids and the correlation of the structures of different classes of materials. In 
the last section of the book the present theoretical problems of luminescence and possible 
future approaches to their solution are discussed. In particular, the need for using 
luminescent solids in the form of single crystals in experimental studies is emphasized. 

This book is not intended to be a comprehensive, up-to-date text-book on luminescence. 
Its value lies in its accounts of the results of recent studies having a rather specialized 
nature and in the discussions which may have more general significance. These discussions 
would probably have been even more valuable if they had been based on abstracts of the 
papers circulated previous to the Conference. It might also have been better to reduce the 
number of papers and to include and expand those in which a marked contribution to the 
knowledge of luminescence processes was evident. 

However, there is included in the book an account of studies showing the renewed 
interest in the release of stored energy in phosphors by infra-red radiation. Quenching and 
stimulation processes in phosphors due to infra-red irradiation were popular subjects for 
early workers in luminescence. Although they are intimately connected with trapped 
electrons, the latest studies show that there is little correlation between the thermal and 
optical ejection of these electrons. The confusion of ideas with respect to the effects of 
infra-red radiation in phosphors shown in the appropriate papers serves to indicate the 
present inadequacy of the theory developed for the interpretation of other processes, such 
as phosphorescence and thermoluminescence. 

Although the book will be of great value to those engaged in research on luminescence, 
the lack of clarity in some of the individual papers will make it difficult reading for others 
having a secondary interest in the subject. G. F. J. GARLICK. 
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ABSTRACTS FOR SECTION B 


An Investigation of the Physical Properties of Thin Films used in the Reduction 
of Surface Reflection, by H. A. DELL. 


ABSTRACT. A simple mathematical treatment of reflection from thin films is extended 
to cover the properties of films likely to be met with in practice. In particular, reflection 
from multible films of more than one substance is considered, and the reflectton of polarized 
light from simple films. A photometer used for testing these results is described, and the 
different types of cryolite films found are enumerated. 

Evidence is given suggesting that polished or aged cryolite films behave as if they 
consisted of a lower stratum of cryolite surmounted by a secondary layer of higher refractive 
index. Some of the properties of this layer are described. 


An Extended Theory of the Reflection Echelon Grating, by L. R. GRIFFIN. 


ABSTRACT. ‘The single reflection echelon grating has the disadvantage that it is 
impossible to interpret the fringe pattern without further information about the correct 
order allocation. The present methods of obtaining this information are briefly reviewed; 
they are progressively less useful as the wavelength decreases into the ultra-violet. The 
properties of the grating are investigated theoretically when the light is incident at relatively 
large angles in two special cases, and the performance of the instrument in these positions 
compared with that of its normal use. It is shown that the required order allocation could 
be obtained using a single reflection echelon grating by measuring the fractional order 
separation with the grating in two positions. This would considerably increase the value 
of the instrument. 


The Duality between Interlinked Electric and Magnetic Circuits and the Formation 
of Transformer Equivalent Circuits, by E. COLIN CHERRY. 


ABSTRACT. When making calculations ona circuit, containing both electric impedances 
and transformers, it is frequently desirable to consider the transformers removed and the 
constraints they impose replaced by a rearrangement of the impedances connected to their 
terminals. Such “equivalent circuits’? may not always be found; the rules are here 
established for their formation, and also for checking, by inspection, whether the transformer 
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constraints are removable in this way in any particular case. It is shown that the 
equivalent electric circuit of a transformer, having any arrangement of magnetic paths, is 
derivable from its magnetic circuit by application of the topological principle of duality. 
This cannot be done if the magnetic circuit is non-planar, as in the case of a transformer 
possessing four or more winding with leakage couplings; a physically realizable circuit does 
not then exist. 

Under certain conditions the principle may be applied in reverse and the impedances in a 
given electric circuit may be coupled together bya suitable transformer, so that the various 
current and voltage constraints are unaltered. 


An Investigation of the Dynamic Elastic Properties of Some High Polymers, by 
K. W. Hitirer and H. Kotsxy. 


ABSTRACT. An apparatus for the investigation of the transmission of sound along 
filaments at frequencies between 1,000 and 6,000 cycles per second is described, both for 
unstrained specimen and whilst they are being elongated at a constant rate of increase 
of strain. Measurements of the dynamic elasticity and damping factors of filaments of 
polythene, neoprene and nylon have been obtained, and the correlation of these results 
with the molecular rearrangements which take place during stretching is discussed. 


A Theory of the Angle of Contact, by G. D. YARNOLD and B. J. Mason. 


ABSTRACT. It does not appear to be possible to give a quantitative account of the angle 
of contact of a continuous liquid in terms of cohesional forces. By using Young’s energy 
relation and making allowance for the effects of adsorbed films on the surface of the solid, 
expressions are derived for the advancing and receding angles, and a general explanation of 
the hysteresis of the angle of contact is put forward. 


The Angle of Contact Between Water and Wax, by G. D. YaRNOLD and 
B. J. Mason. 


ABSTRACT. Theeffects (a) of the velocity of the liquid surface, and (b) of the time of 
immersion, on the angle of contact between clean surfaces of water and wax have been 
investigated by dynamical method. It is clear that Ablett’s measurements are vitiated by 
the fact that the time of immersion is disregarded. The term ‘‘equilibrium angle of contact” 
appears to be meaningless. 


The Absorption of Ultrasonic Waves in Liquids and its Relation to Molecular 
Constitution, by J. M. M. PINKERTON. 


ABSTRACT. Measurements of the absorption of ultrasonic waves in water and ethyl 
alcohol were made by the pulse method, and are more accurate than those made by other 
methods. ~The absorption coefficient a in water was found to vary as the square of the 
frequency v over a wide range of temperatures. ‘The measurements in water were made at 
frequencies between 7:5 and 67:5 Mc/s. and temperatures between 0 and 95°c. The 
absorption decreased with increasing temperature by a factor of 8 between freezing and. 
boiling points. The observed absorption was about three times that calculated from Stokes’ 
formula and this ratio was nearly independent of temperature. In ethyl alcohol measure- 
ments were made at 52:4 Mc/s., and at temperatures between —50° and +60°c. The 
absorption decreased by a factor 3 in this range and observed values were about twice 
those given by Stokes’ formula. Values of (a/v?) at 25° c. are 22-0 X 1071" sec?/cm. for water 
and 50:5 107?’ sec?/cm. for ethyl alcohol. 

Published values of the absorption coefficients of a wide range of pure liquids and of 
some mixtures are reviewed, and a classification proposed. The experimental results. 
suggests that the absorption in associated liquids is due to a different mechanism from that 
in non-associated liquids. In the latter the absorption is probably due to the same pheno- 
menon as in gases: a slow exchange of energy between different degrees of freedom. The 
relaxation frequencies must be above 200 Mc/s. for most liquids. A simple theory is given 
to explain the observed variation of the absorption with concentration in mixtures of non- 
associated liquids. The semi-quantitative argeement of theory with be pt lends. 
support to the relaxation hypothesis. 
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| THE PHYSICAL SOCIETY 


MEMBERSHIP 
Membership of the Society is open to all who are interested in Physics. 


FELLOwsHIP.—A candidate for election to Fellowship must, as a rule, be recommended 
by three Fellows, to two of whom he is known personally. 


StupENT MemsersHip.—A candidate for election to Student Membership must be 
between 18 and 26 years of age and must be recommended from personal knowledge 
by a Fellow. 

Fellows and Student Members may attend all meetings of the Society, and are 
entitled to receive some of the Society’s publications free and others at much reduced 
rates. 

Books and periodicals may be read in the Society’s Library, and a limited number 
of books may be borrowed on application to the Honorary Librarian. 


SUBSCRIPTIONS | 


Fellows pay an Entrance Fee of £1 1s. and an Annual Subscription of £3 3s. 
Student Members pay only an Annual Subscription of 15s. 

Fellows and Student Members may become members of the specialist Groups of 
the Society (Colour, Optical, Low Temperature, Acoustics) without additional sub- 
scription. 


Further information may be obtained from the Secretary-Editor at the Offices of the Society, 
1 LowrHER GARDENS, PRINCE Consort Roap, LONDON S.W. 7 


PHYSICAL SOCIETY SPECIALIST GROUPS 


OPTICAL GROUP 


The Physical Society Optical Group exists to foster interest in and development of all branches o 
optical science. ‘To this end, among other activities, it holds meetings about five times a year to 
discuss subjects covering all aspects of the theory and practice of optics, according to the papers offered. 


COLOUR GROUP 


The Physical Society Colour Group exists to provide an opportunity for the very varied types 0} 
worker engaged on colour problems to meet and to discuss the scientific and technical aspects of their 
work. Five or six meetings for lectures and discussions are normally held each year, and reprints of 
papers are circulated to members when available. A certain amount of committee work is undertaken, 
and reports on Defective Colour Vision (1946) and on Colour Terminology (1948) have already been 
published. 


LOW TEMPERATURE GROUP 


The Low Temperature Group was formed to provide an opportunity for the various groups of 
people concerned with low temperatures—physicists, chemists, engineers, etc.—to meet and become 
familiar with each other’s problems. The Group seeks to encourage investigations in the low temperature 
field and to assist in the correlation and publication of data. 


ACOUSTICS GROUP 

The Acoustics Group was formed to bring together people with varied interests in acoustics—~ 
architectural, musical and biological, as well as physical—so as to advance as rapidly as possible tne 
practice and the theory of the subject. : 


Further information may be obtained from the Offices of the Society : 
1 LowrTHER GARDENS, PRINCE CoNsoRT Roap, LONDON S.W. 7. 


Printed by TayLor AND Francis, Lrp., Red Lion Court, Fleet Street, London E.C.4 


Sa a st i 


